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ABSTRACT 

We  report  on  a  total  of  106  nights  of  optical  interferometric  observations  of  the  e  Aurigae  system  taken  during  the 
last  14  years  by  four  beam  combiners  at  three  different  interferometric  facilities.  This  long  sequence  of  data 
provides  an  ideal  assessment  of  the  system  prior  to,  during,  and  after  the  recent  2009  2011  eclipse.  We  have 
reconstructed  model-independent  images  from  the  10  in-eclipse  epochs  which  show  that  a  disk-like  object  is 
indeed  responsible  for  the  eclipse.  Using  new  three-dimensional,  time-dependent  modeling  software,  we  derive  the 
properties  of  the  F-star  (diameter,  limb  darkening),  determine  previously  unknown  orbital  elements  (D,  /),  and 
access  the  global  structures  of  the  optically  thick  portion  of  the  eclipsing  disk  using  both  geometric  models  and 
approximations  of  astrophysically  relevant  density  distributions.  These  models  may  be  useful  in  future 
hydrodynamical  modeling  of  the  system.  Finally,  we  address  several  outstanding  research  questions  including  mid¬ 
eclipse  brightening,  possible  shrinking  of  the  F-type  primary,  and  any  warps  or  sub-features  within  the  disk. 

Key  words:  binaries:  eclipsing  techniques:  interferometric  techniques:  photometric 
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1.  INTRODUCTION 

Epsilon  (e)  Aurigae  (FKV  0183,  HD  31964)  is  an 
intrinsically  bright  ( V  ~  3.0)  star  system  that  has  confounded 
generations  of  astronomers  since  its  two-year  long,  AV  ^  0.75 
mag  fading  was  first  discovered  nearly  200  years  ago 
(Fritsch  1824).  In  the  early  1900s,  it  was  established  that  e 
Aurigae  was  a  single-line  spectroscopic  binary  with  a  27.1  year 
period  (Ludendorff  1903).  Application  of  the  then  recently 
developed  eclipsing  binary  theory  (Russell  1912a,  1912b)  to  e 
Aurigae  came  to  the  perplexing  result:  the  companion  must  be 
nearly  equal  in  mass  to  the  primary,  but  nearly  invisible  (see 
note  18  of  Shapley  1915).  Although  many  theories  have  been 
proposed  to  explain  this  conundrum  (see  Carroll  et  al.  1991; 
Guinan  et  al.  2002;  Kloppenborg  2012,  for  comprehensive 
summaries  of  the  literature),  it  was  Kopal  (1954)  who  first 
proposed  and  Huang  (1965)  who  later  (independently) 
developed  the  theory  that  the  system  is  a  nearly  edge-on 
eclipsing  binary  composed  of  an  FOIa  supergiant  (primary, 
hereafter  the  F-star)  and  an  unseen  companion  that  was 
believed  to  be  enshrouded  in  a  disk  of  opaque  material.  It  was 
theorized  that  the  passage  of  the  disk  in  front  of  the  F-star 
caused  the  anomalously  long  fading  of  the  system.  This  theory 
received  significant  support  during  the  1984  eclipse  when 
Backman  et  al.  (1984)  detected  the  presence  of  a  500  db  150  K 
blackbody  that  remained  detectable  when  the  F-star  was 
eclipsed. 

The  2009  2011  eclipse  of  e  Aurigae  caused  a  resurgence  of 
interest  in  this  enigmatic  binary,  yielding  a  wealth  of  new 


Max  Planck  Institut  Fur  Radioastronomie. 


scientific  results.  Historic  data  played  a  pivotal  role  in  several 
studies.  Stefanik  et  al.  (2010)  and  Chadima  et  al.  (2010) 
published  two  new,  nearly  identical,  spectroscopic  orbital 
solutions  using  spectroscopic  and  photometric  data  going  back 
to  the  mid-1800s.  Likewise,  Griffin  &  Stencel  (2013) 
conducted  a  meticulous  study  of  historic  and  new  spectroscopic 
observations  covering  the  three  eclipses.  They  discovered  that 
the  precision  by  which  the  spectroscopic  evolution  of  the 
eclipse  unfolds  indicates  that  the  structure  of  the  disk  has  not 
appreciably  changed  on  a  time  scale  of  at  least  the  last  100 
years.  Similarly,  infrared  studies  reported  by  Stencel  et  al. 
(2011)  have  confirmed  persistent  behavior  of  disk  features  over 
at  least  the  last  two  eclipse  cycles. 

New  photometric  and  spectroscopic  observations  have 
revealed  the  progression  of  the  eclipse  in  hitherto  unprece¬ 
dented  detail.  Photometric  monitoring,  largely  by  amateur 
astronomers,  captured  nearly  3700  photometric  points  in  a 
myriad  of  filters  from  the  UV  to  near-IR  (Hopkins  2012). 
Spectroscopic  monitoring  has  shown  that  the  equivalent  width 
of  some  spectral  lines  follow  a  stair  step-like  evolution  pattern 
(Leadbeater  et  al.  2012),  suggesting  some  substructure  exists 
within  the  eclipsing  disk.  For  the  first  time,  a  comprehensive 
spectral  energy  distribution  (SED)  of  the  system  was  assembled 
from  the  UV  to  radio  wavelengths.  These  data  clearly  show  the 
presence  of  not  only  the  F-star  and  a  550  K  disk,  but  also  a 
slight  amount  of  far-UV  flux  thought  to  originate  from  a  B-type 
companion  at  the  center  of  the  disk  (Hoard  et  al.  2010). 
Subsequent  work  has  shown  the  disk  has  an  asymmetric 
temperature  structure  ranging  from  550  1150K  due  to  external 
heating  from  the  F-star  (Hoard  et  al.  2012). 
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In  addition  to  the  classic  observational  methods,  the 
comparatively  new  technique  of  optical  interferometry  was 
recently  applied  to  the  e  Aurigae  eclipse  for  the  first  time.  The 
model-independent  H-band  images  prove,  beyond  any  reason¬ 
able  doubt,  that  an  eclipsing  disk  is  indeed  responsible  for  the 
observed  fadings  (Kloppenborg  et  al.  2010).  Complementary 
spectro-interferometric  observations  of  the  Ha  line  have  shown 
that  the  F-star  has  an  extensive  P-Cygni-like  wind  region  and 
that  the  disk  contains  a  substantial  gaseous  atmosphere  that  can 
eclipse  a  greater  fraction  of  the  F-star  than  the  dark  disk  itself 
(Mourard  et  al.  2012). 

Yet  with  as  much  information  we  have  about  this  system,  we 
have  yet  to  unravel  many  of  its  fundamental  properties. 
Foremost,  the  distance  to  the  system  is  extremely  uncertain, 
ranging  from  0.4  4.0  kpc  as  estimated  by  HIPPARCOS 
(Perryman  &  ESA  1997;  van  Leeuwen  2007).  Dynamical 
parallaxes  from  the  Yerkes  (Strand  1959)  and  Sproul  (van  de 
Kamp  1978;  Heintz  &  Cantor  1994)  observatories  narrow  this 
range  to  0.5  0.7  kpc,  yet  there  is  in  significant  disagreement 
with  the  1.5  ±  0.5  kpc  suggested  via  interstellar  absorption  and 
reddening  (Guinan  et  al.  2012).  Second,  the  evolutionary  state 
of  the  system  has  been  called  into  doubt.  Some  suggest  that  the 
F-star  is  a  massive  (~16  M0)  horizontal  branch  supergiant, 
whereas  others  consider  the  F-star  to  be  a  less  massive  (~2  3 
Mq)  post-AGB  star  (see  Guinan  et  al.  2002;  Hoard  et  al.  2010; 
Sadakane  et  al.  2010,  for  a  review  and  recent  developments 
concerning  this  topic).  Finally,  although  several  models  have 
been  proposed  for  the  eclipsing  disk  (e.g.,  Huang  1965,  1974; 
Wilson  1971;  Takeuti  1986,  2011;  Ferluga  1990;  Lissauer 
et  al.  1996;  Budaj  2011),  these  efforts  were  primarily  based  on 
light  curve  modeling.  Reproducing  the  photometry  is  a 
necessary  condition  of  every  model,  but  is  not  sufficient  for 
proving  the  validity  of  the  hypothesis.  In  particular,  such  work 
cannot  separate  out  the  degeneracies  between  radial  and  height- 
dependent  optical  profiles  in  the  disk. 

We  begin  to  address  several  of  these  issues  by  presenting 
newly  determined  orbital  elements  and  new  models  for  the  disk 
based  upon  simultaneous  photometric  and  interferometric 
modeling  of  the  2009  2011  eclipse  during  the  ingress 
(RJD  =  JD  -  2,400,000  -  55,062  55,193),  totality  (RJD  ~ 
55,193  55,631),  and  egress  (RJD  ~  55,631  55,693)  phases. 
In  Section  2,  we  present  our  observations  and  data  reduction 
methods.  Sections  3  and  4  include  a  discussion  of  our  modeling 
process,  image  reconstruction  methods,  and  statistical  analysis 
of  the  resulting  data.  Results  are  presented  in  Section  5,  and 
finally,  we  draw  conclusions  in  Section  6. 

2.  OBSERVATIONS  AND  DATA  REDUCTION 

This  work  summarizes  a  total  of  106  nights  of  interfero¬ 
metric  observations  taken  by  four  beam  combiners  at  three 
different  interferometric  arrays.  An  account  of  the  observa¬ 
tions  including  the  array,  baselines,  combiner,  spectral 
configuration,  and  calibrators  is  provided  in  Table  1. 
Calibrator  identifiers,  positions,  proper  motions,  uniform 
disk  diameters  (#udd)>  and  the  array  at  which  each  calibrator 
was  used  are  listed  in  Table  2.  The  supplementary  material 
for  this  paper  contain  both  calibrated  and  uncalibrated  data 
saved  in  the  Optical  Interferometry  Exchange  format 
(OIFITS;  Pauls  et  al.  2005).  A  sample  plot  of  the  UV 
coverage,  squared  visibilities,  and  bispectra  for  the  2009-11 
CHARA-MIRC  epoch  is  provided  in  Figure  1.  Equivalent 
plots  for  all  other  epochs  can  be  found  in  the  supplementary 


material,  Figures  Al  A45.  In  the  next  several  paragraphs  we 
describe  the  interferometric  arrays,  combiners,  and  reduction 
methods  in  detail. 

2.7.  Palomar  Testbed  Interferometer  (PTI) 

Our  first  interferometric  data  set  on  e  Aurigae  was  acquired 
using  the  PTI  (Colavita  et  al.  1999)  located  on  Mount  Palomar 
in  California.  The  facility  consisted  of  three  40  cm  siderostats, 
each  located  at  the  termination  of  one  of  the  interferometer’s 
three  arms.  Pairwise  combination  provided  baselines  between 
85  and  110  m.  The  beam  combiner  at  PTI  operated  in  several 
low-resolution  spectral  modes  providing  up  to  11  spectral 
channels  across  the  K  band  (2.2  p m). 

Much  of  the  wide-band  visibility  data  has  been  discussed 
previously  in  Stencel  et  al.  (2008),  therefore,  we  will  not 
consider  it  further.  Here  we  have  re-reduced  the  spectrally 
dispersed  (i.e.,  narrow-band)  data  subject  to  the  calibrator 
diameters  in  Table  2.  We  have  used  the  narrow-  and  wide¬ 
band  calibration  routines  nbCalib  and  wbCalib ,  respectively, 
from  the  V2calib  software  package.  These  programs  are 
available  as  a  web  service,  webCalib ,  on  the  NASA 
Exoplanet  Science  Institute  (formerly  Michelson  Science 
Center)  website.10  We  have  selected  the  PTI  defaults  with  the 
following  exceptions:  (1)  the  calibration  window  was 
extended  to  four  hours;  (2)  no  ratio  correction  was  applied, 
and  (3)  no  minimum  uncertainty  was  enforced  for  reasons 
discussed  in  Stencel  et  al.  (2008).  Output  from  this  pipeline 
were  saved  in  OIFITS  format. 

2.2.  Navy  Precision  Optical  Interferometer  ( NPOI) 

The  NPOI  (Armstrong  et  al.  1998)  is  a  six  telescope  optical 
interferometer  that  started  operation  in  1994.  The  array  may 
be  configured  in  either  an  astrometric  or  imaging  mode.  Data 
from  the  imaging  sub  array  comes  from  six  movable  50  cm 
siderostats  with  baselines  between  16  and  79  m.  The  NPOI 
beam  combiner  operates  at  visible  wavelengths 
(0.5  0.85  pm)  in  16  spectral  channels.  The  NPOI  observa¬ 
tional  setup  and  data  recording  procedure  can  be  found  in 
Hummel  et  al.  (2003)  and  Benson  et  al.  (2003).  Post¬ 
processing  and  data  reduction  were  performed  using  C. 
Hummel’s  OYSTER  software  package. 

The  NPOI  observed  e  Aurigae  on  a  total  of  29  nights 
between  2006  February  and  2010  April  as  shown  in  Table  1. 
The  data  were  initially  calibrated  with  respect  to  HD  32630 
assuming  a  uniform  disk  diameter  of  0.507  ±  0.025 
milliarcseconds  (mas)  and  saved  as  OIFITS  files.  Prior  to 
modeling  these  data,  we  recalibrated  the  OIFITS  files  by 
multiplying  the  visibilities  and  closure  amplitudes  by  the  ratio 
of  the  uniform  disk  function  of  the  former  and  new  diameter  as 
listed  in  Table  2. 

2.3.  Center  for  High  Angular  Resolution  Astronomy  ( CHARA ) 

Georgia  State  University’s  CHARA  (ten  Brummelaar 
et  al.  2005)  is  an  interferometric  array  located  on  Mount 
Wilson,  CA.  The  array  consists  of  six  1-m  telescopes  that  can 
be  combined  to  form  up  to  15  baselines  ranging  in  length  from 
34  to  331  m.  Using  the  longest  baselines,  a  resolution  of  down 
to  0.5  mas  in  the  H  band  (0.7  mas  in  the  K  band)  can  be 
realized. 


10  http://nexsciweb.ipac.caltech.edu/cgi  bin/webCalib/webCalib.cgi 
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Table  1 

List  of  Observations  from  the  MIRC,  NPOI,  CLIMB,  and  PTI  Beam  Combiners 


Date 

JD 

Configuration 

Array 

Combiner /Mode 

Calibrators  (HD) 

Comment 

1997  10  22 

2,450,744.02 

N  S 

PTI 

K  band,  7  channels 

32630,  33167 

1997  11  09 

2,450,761.95 

N  S 

PTI 

K  band,  1 1  channels 

32630 

1998  11  07 

2,451,124.95 

N  S 

PTI 

K  band,  5  channels 

32630 

1998  11  25 

2,451,142.94 

N  S 

PTI 

K  band,  5  channels 

30823 

1998  11  26 

2,451,143.93 

N  S 

PTI 

K  band,  5  channels 

42807 

2005  12  11 

2,453,715.78 

N  W 

PTI 

K  band,  9  channels 

29645 

2006  01  31 

2006  02  25 

2,453,791.74 

E06  AC0  AW0, 

PTI 

NPOI 

Kband 

16  channels 

no  calibrators 

32630 

2006  02  26 

2,453,792.75 

AC0  AE0  AN0 

E06  AC0  AW0, 

NPOI 

16  channels 

32630 

2006  03  03 

2,453,432.68 

AC0  AE0  AN0 

E06  ACO  AWO, 

NPOI 

16  channels 

32630 

Noisy  data  on  longest  baseline 

2007  03  02 

2,454,161.63 

ACO  AEO  ANO 

W07  ACO  ANO, 

NPOI 

16  channels 

32630 

2007  03  05 

2,454,164.67 

AWO  ACO  W07 

E06  ANO  AWO, 

NPOI 

16  channels 

32630 

2007  03  10 

2,454,169.65 

E06  AWO  W07 

E06  ANO  AWO, 

NPOI 

16  channels 

32630 

2007  03  12 

2,454,898.63 

E06  AWO  W07 

NPOI 

16  channels 

32630 

Calibrator  scans  incoherent,  unusable 

2007  03  13 

2,454,172.63 

E06  ANO  AWO, 

NPOI 

16  channels 

32630 

2007  03  14 

2,454,173.63 

E06  AWO  W07 

E06  ANO  AWO, 

NPOI 

16  channels 

32630 

2007  03  15 

2,454,174.63 

E06  AWO  W07 

E06  ANO  AWO, 

NPOI 

16  channels 

32630 

2007  03  16 

2,454,175.64 

E06  AWO  W07 

E06  ANO  AWO, 

NPOI 

16  channels 

32630 

2007  10  19 

2,454,392.99 

E06  AWO  W07 

N  S 

PTI 

K  band,  5  channels 

29645 

2007  10  20 

2,454,393.97 

N  S 

PTI 

K  band,  9  channels 

30138,  32630 

Incoherent  v2 

2007  10  21 

2,454,394.98 

N  S 

PTI 

K  band,  5  channels 

30138 

2007  11  27 

2,454,431.85 

N  S 

PTI 

K  band,  5  channels 

27946,  30138,  32537 

2007  12  23 

2,454,457.78 

N  W 

PTI 

K  band,  5  channels 

30138,  32630,  33167 

Incoherent  v2 

2007  12  24 

2,454,458.75 

N  W 

PTI 

K  band,  5  channels 

32630 

2008  02  16 

2,454,512.63 

N  W 

PTI 

K  band,  5  channels 

30138,  32630 

2008  02  17 

2,454,513.64 

N  W 

PTI 

K  band,  5  channels 

30138,  32630 

2008  02  18 

2,454,514.65 

N  W 

PTI 

K  band,  5  channels 

30138,  32630 

2008  09  19 

2,454,729.02 

SI  El  W1  W2 

CHARA 

MIRC,  4T,  LR  H, 

3360,  22928,  219080 

2008  10  17 

2,454,757.00 

N  W 

PTI 

choppers 

K  band,  5  channels 

32630 

2008  10  26 

2,454,765.98 

N  S,  N  W 

PTI 

K  band,  5  channels 

29645,  30138,  30823,  32630 

Incoherent  v2 

2008  11  07 

2,454,777.89 

SI  El  W1  W2 

CHARA 

MIRC,  4T,  LR  H, 

5448,  101501,  24398,  50019 

2008  11  08 

2,454,778.88 

SI  El  W1  W2 

CHARA 

choppers 

MIRC,  4T,  LR  H, 

5448,  101501,  24398,  50019 

2008  11  08 

2,454,778.87 

N  S 

PTI 

choppers 

K  band,  5  channels 

29645,  30138,  30823,  32630 

2008  11  09 

2,454,779.87 

N  W 

PTI 

K  band,  5  channels 

29645,  30138,  30823,  32630 

2008  11  15 

2,454,786.87 

N  S 

PTI 

K  band,  5  channels 

42807,  73262 
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Date 

JD 

Configuration 

Array 

Combiner /Mode 

2008  11  22 

2,454,792.82 

N  S,  N  W 

PTI 

K  band,  5  channels 

2008  12  10 

2,454,810.82 

SI  El  W1  W2 

CHARA 

MIRC,  4T,  LR  H, 

choppers 

2008  12  20 

PTI 

2008  12  21 

PTI 

2009  03  08 

2,454,898.63 

E06  ANO  AWO, 

E06  AWO  W07 

NPOI 

16  channels 

2009  03  12 

2,454,902.63 

E06  ANO  AWO, 

E06  AWO  W07 

NPOI 

16  channels 

2009  11  02 

2,455,137.80 

SI  El  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

2009  11  02 

2,455,137.95 

S2  E2  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

2009  11  03 

2,455,138.79 

SI  El  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

2009  11  03 

2,455,138.98 

S2  E2  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

2009  11  04 

2,455,139.75 

SI  El  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

2009  11  04 

2,455,139.93 

S2  E2  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

2009  12  02 

2,455,167.77 

SI  El  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

2009  12  02 

2,455,167.89 

S2  E2  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

2009  12  03 

2,455,168.74 

SI  El  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

2009  12  03 

2,455,168.93 

S2  E2  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

2009  12  04 

2,455,169.88 

S2  E2  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

2009  12  18 

2,455,183.77 

AWO  ACO,  AEO  ACO 

NPOI 

16  channels 

2009  12  19 

2,455,184.73 

AWO  ACO,  AEO  ACO 

NPOI 

16  channels 

2009  12  20 

2,455,185.80 

AWO  ACO,  AEO  ACO 

NPOI 

16  channels 

2009  12  21 

2,455,186.81 

ACO  E06  AWO, 

AEO  ACO  ANO 

NPOI 

16  channels 

2009  12  25 

2,455,190.71 

AWO  ACO,  AEO  ACO 

NPOI 

16  channels 

2009  12  27 

2,455,192.69 

AWO  ACO,  AEO  ACO 

NPOI 

16  channels 

2010  01  03 

2,455,199.67 

ACO  E06  AWO, 

AEO  ACO  ANO 

NPOI 

16  channels 

2010  01  04 

2,455,289.60 

NPOI 

16  channels 

2010  01  05 

2,455,201.68 

ACO  E06  AWO, 

AEO  ACO  ANO 

NPOI 

16  channels 

2010  01  06 

2,455,202.79 

AEO  ANO  AWO, 

AWO  AEO  ANO 

NPOI 

16  channels 

2010  01  07 

2,455,203.78 

AEO  ANO  AWO, 

AWO  AEO  ANO 

NPOI 

16  channels 

2010  01  08 

2,455,204.66 

AEO  ANO  AWO, 

AWO  AEO  ANO 

NPOI 

16  channels 

2010  01  09 

2,455,205.86 

AEO  ANO  AWO, 

AWO  AEO  ANO 

NPOI 

16  channels 

2010  01  10 

2,455,206.83 

AEO  ANO  AWO, 

AWO  AEO  ANO 

NPOI 

16  channels 

2010  01  12 

2,455,208.77 

AEO  ANO  AWO, 

NPOI 

16  channels 

AWO  AEO  ANO 
AEO  ANO  AWO, 
AWO  AEO  ANO 


2010  01  13 


2,455,209.70 


NPOI 


16  channels 


32630 
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Table  1 

(Continued) 


Date 

JD 

Configuration 

Array 

Combiner/Mode 

Calibrators  (HD) 

Comment 

2010  01  16 

2,455,212.69 

AC0  E06  AW0, 

AE0  AC0  AN0 

NPOI 

16  channels 

32630 

2010  02  14 

2,455,241.62 

AC0  E06  AW0, 

AE0  ACO  ANO 

NPOI 

16  channels 

32630 

2010  02  15 

2,455,242.62 

ACO  E06  AWO, 

AEO  ACO  ANO 

NPOI 

16  channels 

32630 

2010  02  16 

2,455,243.62 

ACO  E06  AWO, 

AEO  ACO  ANO 

NPOI 

16  channels 

32630 

2010  02  17 

2,455,244.62 

ACO  E06  AWO, 

AEO  ACO  ANO 

NPOI 

16  channels 

32630 

2010  02  18 

2,455,245.74 

S2  E2  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

32630,  41636 

2010  02  19 

2,455,246.70 

SI  El  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

No  calibrators,  used  closure 

See  Text 

phase  only 

2010  04  03 

2,455,289.60 

ANO  E06  AWO, 

AWO  E06  W07 

NPOI 

16  channels 

32630 

2010  08  20 

2,455,428.96 

SI  El  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

3360,  21770,  24760,  32630 

2010  08  21 

2,455,429.96 

S2  E2  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

3360,  21770,  24760,  32630 

2010  08  22 

2,455,430.99 

SI  El  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

21770,  32630 

2010  08  23 

2,455,431.96 

S2  E2  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

3360,  21770,  24760,  32630, 

219080 

2010  09  23 

2,455,462.92 

SI  E2  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

24760,  32630 

2010  09  24 

2,455,463.92 

S2  E2  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

3360, 6961,  21770,  219080,  24760, 

32630 

2010  09  26 

2,455,465.98 

W1  W2  S2  E2 

CHARA 

MIRC,  4T,  LR  H 

32630 

MIRC  +  PAVO 

2010  09  27 

2,455,466.99 

W1  W2  S2  E2 

CHARA 

MIRC,  4T,  LR  H 

32630 

MIRC  +  PAVO 

2010  09  28 

2,455,467.97 

W1  W2  E2  S2 

CHARA 

MIRC,  4T,  LR  H 

32630 

MIRC  +  PAVO 

2010  10  26 

2,455,495.85 

SI  El  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

24760,  32630 

2010  10  26 

2,455,495.98 

S2  E2  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

32630 

2010  10  27 

2,455,496.86 

SI  El  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

32630 

2010  10  27 

2,455,496.91 

S2  E2  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

32630 

2010  11  04 

2,455,504.95 

S2  E2  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

32630 

2010  11  05 

2,455,505.89 

SI  El  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

21770,  32630,  41636 

2010  12  12 

2,455,542.77 

W1  W2  E2  SI 

CHARA 

MIRC,  4T,  LR  H 

32630,  50019 

MIRC  +  PAVO 

2010  12  13 

2,455,543.75 

W1  W2  E2  SI 

CHARA 

MIRC,  4T,  LR  H 

24760,  32630 

MIRC  +  PAVO 

2010  12  14 

2,455,544.75 

W1  W2  S2  El 

CHARA 

MIRC,  4T,  LR  H 

32630 

MIRC  +  PAVO 

2011  01  18 

2,455,579.67 

SI  El  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

32630 

2011  01  18 

2,455,579.89 

S2  E2  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

32630,  50019 

2011  01  19 

2,455,580.66 

SI  El  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

32630 

2011  01  19 

2,455,580.78 

S2  E2  W1  W2 

CHARA 

MIRC,  4T,  LR  H 

32630,  41636 

2011  03  18 

El  E2  W1 

CHARA 

CLIMB,  H 

32630 

Non  standard  readout  mode  and  bad  conditions.  Not  usable. 

2011  04  01 

2,455,652.68 

El  E2  W1 

CHARA 

CLIMB,  K 

32630 

2011  04  03 

2,455,654.67 

S2  W2  W1 

CHARA 

CLIMB,  K 

32630 

2011  04  05 

2,455,656.66 

S2  El  W1 

CHARA 

CLIMB,  K 

32630 

2011  09  18 

2,455,823.03 

W1  S2  SI  El  E2  W2 

CHARA 

MIRC,  6T,  LR  H 

27396 

2011  09  24 

2,455,829.03 

W1  S2  SI  El  E2  W2 

CHARA 

MIRC,  6T,  LR  H 

32630 

See  Text 

2011  10  10 

2,455,844.94 

W1  S2  SI  El  E2  W2 

CHARA 

MIRC,  6T,  LR  H 

32630 

2011  11  03 

2,455,868.85 

W1  S2  SI  El  E2  W2 

CHARA 

MIRC,  6T,  LR  H 

21770,  24760,  32630 

Note.  Calibrator  IDs  may  be  cross  referenced  with  Table  2. 
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Table  2 

Calibrators  and  Adopted  Uniform  Disk  Diameters  (#udd) 


HD 

Name 

Position  (J2000) 

(mas  yr  l) 

Ms 

(mas  yr  ) 

7 r 

(mas) 

Published  diameters 

Adopted 

Array 

References  Notes 

R.A. 

(HH  MM  SS.SS) 

decl. 

(DD  MM  SS.SS) 

#UDD-H 

(mas) 

#UDD-K 

(mas) 

°"0UDD 

(mas) 

#UDD 

(mas) 

^UDD 

(mas) 

3360 

zet  Cas 

00  36  58.28419 

+53  53  48.8673 

17.38 

9.86 

5.5 

0.287 

0.288 

0.02 

0.287 

0.020 

C 

(2) 

5448 

37  And 

00  56  45.21211 

+38  29  57.6380 

153.48 

36.49 

25.14 

0.593 

0.594 

0.042 

0.593 

0.042 

C 

(2) 

6961 

tet  Cas 

01  11  06.16225 

+55  08  59.6472 

226.77 

18.75 

24.42 

0.471 

0.033 

c 

(4) 

21770 

36  Per 

03  32  26.26028 

+46  03  24.6965 

52.11 

75.26 

27.53 

0.582 

0.041 

c 

(4) 

22928 

del  Per 

03  42  55.50426 

+47  47  15.1746 

25.58 

43.06 

6.32 

0.549 

0.038 

c 

(4) 

24398 

zet  Per 

03  54  07.92248 

+31  53  01.0812 

5.77 

9.92 

4.34 

0.700 

0.030 

c 

(1) 

27396 

53  Per 

04  21  33.16557 

+46  29  55.9554 

20.06 

35.45 

6.43 

0.285 

0.020 

c 

(2) 

27946 

67  Tau 

04  25  25.01518 

+22  11  59.9876 

111.97 

47.71 

22.03 

0.442 

0.443 

0.031 

0.443 

0.031 

p 

(2) 

28052 

71  Tau 

04  26  20.74092 

+15  37  05.7652 

114.31 

32.19 

20.37 

0.542 

0.038 

p 

(4) 

29645 

HR  1489 

04  41  50.25660 

+38  16  48.6622 

241.65 

97.15 

31.38 

0.521 

0.522 

0.037 

0.523 

0.037 

p 

(2) 

30138 

HR  1514 

04  46  44.47871 

+40  18  45.3270 

9.07 

36.83 

7.53 

0.856 

0.86 

0.061 

0.826 

0.061 

p 

(4) 

30823 

HR  1550 

04  52  47.75706 

+42  35  11.8569 

10.63 

0.46 

7.2 

0.317 

0.317 

0.022 

0.317 

0.022 

p 

(2)  Get  cal  estimates  ~0.1  mas  larger 

32537 

9  Aur 

05  06  40.62967 

+51  35  51.8025 

30.49 

172.89 

38.04 

0.589 

0.041 

p 

(4) 

32630 

eta  Aur 

05  06  30.89337 

+41  14  04.1127 

31.45 

67.87 

13.4 

0.453 

0.012 

C,P,N 

(3) 

33167 

HR  1668 

05  10  42.92081 

+46  57  43.4550 

58.34 

149.99 

20.54 

0.498 

0.035 

P 

(4) 

37147 

122  Tau 

05  37  03.73543 

+17  02  25.1776 

42.24 

33.69 

20.58 

0.375 

0.026 

P 

(4) 

41636 

HR  2153 

06  08  23.13611 

+41  03  20.6194 

1.36 

48.55 

7.82 

0.765 

0.054 

c 

(4) 

42807 

HR  2208 

06  13  12.50242 

+10  37  37.7095 

77.38 

298 

55.71 

0.473 

0.475 

0.034 

0.475 

0.034 

p 

(2) 

50019 

tet  Gem 

06  52  47.33887 

+33  57  40.5175 

1.66 

47.31 

17.25 

0.802 

0.804 

0.056 

0.802 

0.056 

c 

(2) 

73262 

del  Hya 

08  37  39.36627 

+05  42  13.6057 

70.19 

7.9 

20.34 

0.464 

0.465 

0.032 

0.465 

0.032 

p 

(2) 

101501 

61  UMa 

11  41  03.01636 

+34  12  05.8843 

12.55 

380.75 

104.04 

0.864 

0.061 

c 

(4) 

219080 

7  And 

23  12  33.00380 

+49  24  22.3455 

90.23 

95.56 

40.67 

0.665 

0.047 

c 

(4) 

Notes.  C:  CHARA,  N:  NPOI,  P:  PTI. 

References  (1)  MIRC  Calibrator  database,  (2)  Lafrasse  et  al.  (2010),  (3)  Maestro  et  al.  (2013),  (4)  Value  computed  using  SearchCal. 
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Initial  calibration  observations  of  e  Aurigae  were  taken  far  in 
advance  of  the  eclipse,  in  2008  October/November/December. 
Semi-regular  observations  were  scheduled  around  the  photo¬ 
metric  eclipse,  beginning  in  2009  October /November  and 
ending  in  2011  November.  In  total,  e  Aurigae  was  observed  on 
38  nights,  yielding  19  individual  epochs  after  consecutive 
nights  were  merged.  The  first  two  eclipse  ingress  epochs  were 
previously  discussed  in  Kloppenborg  et  al.  (2010)  in  which 
reconstructed  images  and  a  preliminary  model  (comprised  of  an 
infinitely  thin,  but  optically  thick  disk  seen  in  projection)  for 
the  eclipsing  disk  were  presented. 

2.3.1.  Michigan  InfraRed  Combiner  ( MIRC ) 

MIRC  (described  in  Monnier  et  al.  2004;  Monnier  2006) 
went  through  several  revisions  during  our  observing  program. 
As  it  was  first  used,  the  combiner  was  configured  for  four- 
telescope  beam  combination  using  a  camera  sensitive  in  the 
H  and  K  bands  with  three  spectral  resolution  options  (R  ~  44, 
150,  or  400).  We  followed  the  standard  observing  procedures 
and  reduced  the  data  using  the  data  reduction  pipeline 
described  in  Monnier  et  al.  (2007).  After  co-adding  frames, 
background  subtraction,  and  a  Fourier  transform  of  the  raw 
data,  fringe  amplitudes  and  phases  are  used  to  form  the  squared 
visibilities  and  triple  products.  These  data  are  written  to 
OIFITS  files  for  further  analysis. 

In  the  earliest  data,  photometric  calibration  was  achieved 
using  choppers  that  temporally  encoded  the  flux  coming  from 
each  telescope  by  periodically  blocking  each  beam  at  a  unique 
frequency.  In  2009  August,  the  choppers  were  replaced  by 
dedicated  photometric  channels  (Che  et  al.  2010)  that  utilize  a 
fraction  of  light  from  each  science  beam  for  calibrated, 
spectrally  dispersed  photometry.  This  led  to  a  dramatic 
improvement  in  uncertainties.  In  the  spring  of  2011,  MIRC 
was  upgraded  again  to  combine  light  from  all  six  of  CHARA’s 
telescopes  (see  initial  report  in  Monnier  et  al.  2010),  permitting 
measurement  of  15  non-redundant  visibilities  and  all  10 
independent  closure  phases  available  at  CHARA. 

2.3.2.  CLassic  Interferometry  with  Multiple  Baselines  (CLIMB) 

During  the  egress  phase  of  the  2009  2011  eclipse,  when 
MIRC  was  unavailable  for  use,  we  employed  the  CLIMB 
(Sturmann  et  al.  2010)  combiner.  CLIMB  is  a  three-telescope 
beam  combiner  which  operates  in  one  of  five  broadband 
spectral  modes.  For  our  observations  we  configured  CLIMB  in 
with  the  K  band  (A  =  2.133  //m,  FWHM  0.349 /xm)  filter 
installed  and  destructive  readout  mode.  Because  of  the  large 
hour  angle  during  our  observations,  the  UV  coverage  afforded 
by  these  data  is  quite  limited. 

We  used  the  standard  observing  technique  and  reduced  our 
data  with  the  CLIMB  version  2.1  reduction  pipeline  (ten 
Brummelaar  et  al.  2012).  Data  were  calibrated  subject  to  the 
calibrators  listed  in  Table  2  and  written  to  OIFITS  format  for 
further  analysis. 

3.  MODELING  AND  ANALYSIS 

When  our  interferometric  observations  at  CHARA  began,  no 
publicly  available  software  was  capable  of  modeling  both 
interferometric  and  photometric  data  in  a  fully  3D,  time 
dependent  fashion.  We  wrote  two  software  packages  to  solve 
this  problem:  The  Simulation  and  Modeling  Tool  for  Optical 
Interferometry  (SIMTOI)  and  the  OpenCL  Interferometry 


Library  (liboi).  Because  there  are  no  prior  publications 
discussing  these  software,  we  present  a  brief  overview  of  their 
capabilities  here. 

SIMTOI 11  (Kloppenborg  &  Baron  2012b)  is  an  open  source 
(GPL)  C/C++  program  for  fitting  time-dependent,  3D  models 
to  large  data  sets.  Instead  of  using  analytic  models,  SIMTOI 
uses  the  Open  Graphics  Library  (OpenGL)  to  tessellate  the 
“surfaces”  of  the  astronomical  objects  being  simulated. 
Geometrical  effects,  such  as  limb  darkening,  are  applied  using 
programs  written  in  the  OpenGL  Shading  Language  (GLSL). 
The  models  may  be  positioned  statically  in  their  3D 
environment  using  fixed  Cartesian  coordinates  or  be  deter¬ 
mined  dynamically  from  a  Keplerian  orbit.  All  models  also 
have  an  intrinsic  inclination,  i,  position  angle,  Q,  and  rotational 
zero  point,  uj.  The  geometrical  primitives  are  rotated,  translated, 
and  orthographically  rendered  to  a  multisample  anti-aliasing 
buffer.  OpenGL  computes  the  total  flux  in  each  pixel  using  a 
method  akin  to  ray  tracing,  creating  2D  images  against  which 
observed  data  is  compared.  At  present,  SIMTOI  supports  both 
interferometric  and  photometric  data.  Additional  data  types  can 
be  supported  by  subclassing  and  registering  the  new  data  type 
with  SIMTOT  s  plugin  API.  SIMTOI  provides  a  similar  API  for 
registering  minimization  engines.  Presently  implemented 
minimizers  include  a  recursive  grid  search,  Levmar  (Lour- 
akis  2005),  MultiNest  (Feroz  &  Hobson  2008;  Feroz 
et  al.  2009),  and  a  bootstrapping  minimizer  based  on  Levmar. 
The  validity  of  models  produced  by  SIMTOI  has  been  tested 
against  iota  Peg,  the  CHARA-MIRC  closure  phase  calibrator 
(Monnier  et  al.  2007)  and  LitPro  (Tallon-Bosc  et  al.  2008) 
model  results. 

SIMTOI  uses  liboi 12  (Kloppenborg  &  Baron  2012a)  to 
generate  interferometric  observables,  liboi  is  an  open-source 
(LGPL)  C/C++  library  which  implements  the  backend  from 
the  GPu  Accelerated  Image  Reconstruction  (Baron  &  Klop¬ 
penborg  2010)  program,  liboi  aims  to  provide  software 
developers  with  convenient  access  to  fast  routines  for  common 
interferometric  tasks.  The  software  heavily  relies  on  the 
heterogeneous  computing  environment  of  the  Open  Compute 
Language  (OpenCL)  to  target  a  wide  range  of  traditional  and 
multi-core  CPUs;  servers,  hand-held/embedded  devices,  spe¬ 
cialized  hardware,  and  Graphical  Processing  Units  (GPUs). 
The  computational  correctness  of  liboi  is  provided  by  a  series 
of  built-in  unit  tests  to  analytical  functions. 

3.1.  Bayesian  Model  Selection 

In  this  work,  we  use  Bayesian  statistics  to  assess  the  relative 
goodness  of  fit  between  our  proposed  models  rather  than 
traditional  chi-squared  methods  that  are  ill-adapted  to  make 
such  inferences  (Marshall  et  al.  2006).  Bayesian  statistics 
provides  a  consistent  approach  to  estimate  a  set  of  parameters, 
0,  in  a  hypothesis  (e.g.,  an  image  or  model),  H ,  given  some 
observed  data,  D.  Bayes’  theorem  states  that 

f(e|P.«)=f<°'e-">f(e'">.  a) 

P(D\H ) 

where  P(0\D,  H)  =  P(Q)  is  the  posterior  probability 
distribution  of  the  parameters,  P(D\Q,  H)  =  L(0)  is  the 
likelihood,  P(Q\H )  =  7r(0)  is  the  prior,  and  P(D\H)  =  Z  is 
the  Bayesian  Evidence.  In  parameter  estimation  problems, 

1 1  SIMTOI:  https: / / github.com/bkloppenborg / simtoi 

12  liboi :  https://github.com/bkloppenborg/liboi 
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Figure  1.  UV  coverage,  data,  beat  fit  bootstrapped  model,  and  residual  plot  for  the  2009  1 1  CHARA  MIRC  epoch.  These  data,  taken  during  the  ingress  phase  of  the 
eclipse,  shew  a  clear  departure  from  circular  symmetry  in  all  measured  quantities  as  evidenced  by  the  visibilities,  triple  amplitudes,  and  closure  phases  being 
significantly  different  at  similar  baseline/triplet  lengths  (baseline  lengths  summed  in  quadrature).  Equivalent  figures  for  other  epochs  can  be  found  in  extended  Hgures 
A1  A45. 

(The  complete  figure  set  (45  images)  is  available.) 


where  the  model  remains  the  same,  the  normalization  factor,  Z, 
is  often  ignored  as  it  is  independent  of  the  parameters  0.  When 
selecting  between  various  models,  the  evidence  plays  a  central 
role  through  the  Bayes  Factor: 


A  principal  difficulty  in  using  Bayesian  evidence  for  model 
selection  is  that  the  multidimensional  integral, 

P(D\H)  =  Z=  jL(Q)n(©)dD©,  (3) 


j;  l  fol”)  ^  0 

/>(/70|D)  P(D\Ho)P{H0)  Z0P{Ho)' 

The  evidence  is  the  average  of  the  likelihood  over  the  prion 
Thus  a  simple  model  with  greater  likelihood  over  the  parameter 
range  will  be  favored  over  a  more  complex  model  with  lower 
likelihood  over  the  parameter  range,  unless  the  latter  is 
significantly  better  at  explaining  the  data.  Therefore  the  Bayes 
factor  automatically  implements  the  Occam  razor  principle. 


must  be  evaluated.  For  this  work,  we  decided  to  use  the 
MultiNest  library  (Feroz  &  Hobson  2008;  Feroz 
et  al.  2009,  2013)  that  numerically  estimates  this  integral  by 
intelligently  exploring  the  parameter  range  using  Markov  chain 
methods  and  ellipsoidal  bounding  conditions.  In  the  model 
fitting  process  we  have  assumed  non-inf ormative  (flat)  priors 
and  used  the  standard  likelihood  function: 


£(©)  =  n  fd—  exp 

t 


(•V/ ! -  Djf 
2a? 


(4) 
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where  M{9)t  is  the  model  prediction,  D,  is  the  data,  and  pf  is  the 
uncertainty  associated  for  the  /th  data  point, 

3.2.  F-star  Model 

Within  the  SIMTOI  framework,  the  F-star  is  modeled  as  a 
single  uniformly  illuminated  sphere,  located  at  the  origin,  to 
which  GLSL  shaders  (simulating  limb  darkening)  were  applied. 
The  PT1  observations  are  very  high  up  on  the  visibility  curve, 
therefore,  they  provide  no  reliable  measurement  of  limb 
darkening.  Consequently,  we  fit  all  PT1  data  with  a  uniform 
disk  model.  The  CHARA  and  NPOl  observations  frequently 
resolve  the  F-star  near  or  beyond  the  first  visibility  null.  These 
data  show  clear  departures  from  uniform  disk  behavior.  To 
account  for  this,  we  implemented  several  limb  darkening  laws 
in  GLSL.  These  include  linear,  logarithmic,  square  root,  power 
law  (via.  Hestroffer  1997),  and  a  few  multi-parameter  laws 
(e.g..  Claret  2000;  Claret  Si  Hauschildt  2003;  Fields 
et  al.  2003). 

3.3.  Disk  Mod t1  Is 

The  previous  models  for  the  c.  Aurigae  disk  were  purely 
geometric  representations  with  hard  edges  (e.g.,  Kloppenborg 
et  al,  2010;  Kloppenborg  2012,  and  references  therein).  In  this 
work,  we  have  created  geometric  and  astrophysical  density 
distribution  models  with  position-dependent  optical  properties. 
Because  of  the  edge-on  nature  of  the  eclipse,  we  elected  to 
represent  all  disks  as  a  series  of  concentric  ri  ngs  of  infinitesimal 
thickness  and  uniform  total  height,  h ,  that  are  equally  spaced 
between  an  inner.  rm7  and  outer,  rou{,  radius.  The  rings  are 
connected  at  the  midplane  by  another  surface. 

Inspired  by  the  appearance  of  proplyds  in  the  Orion  Nebula 
(e.g.,  Ricci  et  al,  2008),  the  geometrical  models  have  opacity 
(via.  OpenGL  source  transparency,  srca)  that  is  controlled  by  a 
double  power  law  that  is  a  function  of  radius  and  height: 


to  be  ascertained,  we  have  treated  the  product  kqPq  as  a  single 
quantity  in  our  minimization. 

In  this  work  we  also  introduce  time -dependence  by 
positioning  the  disk  according  to  a  Keplerian  orbit.  The 
spectroscopic  orbital  solutions  from  Stefan ik  et  al.  (2010)  and 
Chadima  et  al.  (2010)  are  in  excellent  agreement  with  each 
other;  however,  the  astrometric  orbital  parameters  are  not  well 
constrained.  Therefore,  our  model  adopts  the  spectroscopic 
orbital  parameters  (longitude  of  periastron,  w,  eccentricity,  e, 
period,  F,  and  time  of  periastron,  T)  from  Stefan  ik  et  al,  (2010) 
and  we  derive  the  astrometric  parameters  (inclination,  i, 
longitude  of  the  ascending  node,  fi,  and  total  orbital  semimajor 
axis,  Qj)  via.  minimization  to  our  data. 

3.4.  Modeling  Process 

Within  the  framework  of  the  aforementioned  models,  the 
number  of  parameters  needed  ranges  from  6  in  the  case  of  a 
uniformly  illuminated  F-star  and  cylindrical  disk  (e.g.,  i,  Q,  a?, 
$itdd.  'disk.  and  frdisk)  to  10  with  a  limb  darkened  disk  and  more 
complex  opacity  model  for  the  disk  (see  Table  4  for  a  summary 
of  all  parameters  used  and  their  permitted  ranges).  To  reduce 
the  number  of  degrees  of  freedom,  we  elected  to  establish 
bounds  by  solving  a  series  of  subproblems  first,  then  lift 
restrictions  to  generalize  the  results. 

1.  First,  use  Bayesian  evidence  and  the  post -eclipse 
CHARA-M1RC  data  to  establish  bounds  for  the  diameter 
and  limb-darkening  of  the  F-star, 

2,  Once  determined,  the  diameter  was  used  to  approximate 
the  total  orbital  semi  major  axis,  aj,  for  the  system  using 
two  methods.  To  first  order,  we  may  assume  the  orbit  is 
circular  and  the  eclipse  transects  the  equator  of  the  F-star. 
Then  oj  may  be  found,  to  first  order,  by  equating  the 
fraction  of  the  orbit  spent  in  ingress,  with  the  equivalent 
sector  of  the  orbit  via. 


fe)  ■  <5> 

For  the  two  astrophysical  disk  models,  we  keep  the 
concentric  ring  representation  of  the  disk,  but  modify  the 
opacity  of  each  ring  according  to  a  real  density  distribution. 
The  first  model  from  Pascucci  et  al.  (2004)  consists  of  a  power 
law  in  the  radial  direction  and  a  scale  height  exponential  taper 
in  the  vertical  direction: 


=  (8) 
^ingress  s  *4tar 

where  p  is  the  perimeter  of  the  orbit,  s  is  the  orbital 
sector,  aT  is  the  separation  of  the  components,  and  #star  is 
the  diameter  of  the  F-star  in  radians.  If  s  <£  p,  we  may 
perform  a  second -order  approximation  for  the  sector 
using  the  Ramanujan  approximation  for  the  perimeter  of 
an  ellipse.  In  this  case,  we  find: 


pir,  ^)  =  j  e  dh)  ,  (6) 

where  p0  is  the  density,  rc  is  the  scale  radius,  h  =  hc  hc 

is  the  scale  height,  and  units  are  in  angular  quantities  (e.g.,  mas 
or  mas  3  when  appropriate).  The  second  density  distribution  is 
characterized  by  a  power-law  in  the  inner  disk  and  an 
exponential  taper  at  large  radii  (cf.  Andrews  et  al.  2009): 


z)  =  Po 


1 

2 


(7) 


For  these  models,  the  OpenGL  transparency  is  calculated  as 
srca  =  1  —  to  ensure  the  pixel  intensity  is  com¬ 

puted  following  radiative  transfer  conventions.  Because  the 
composition  and  precise  optical  properties  of  the  disk  have  yet 


— - —  ra  ttcljT 3 f  1  +  V 1  ~ 

ingress 

-  JloVl  -e2  +3(2  -*2)  1^4  (9) 

i  Pstar 

where  uj  and  e  are  the  aforementioned  orbital  quantities. 

3.  Approximate  the  position  angle  of  the  ascending  node,  fi, 
by  computing  the  average  of  the  position  angles 
determined  from  single-epoch  in-eclipse  minimizations, 

4.  Determine  the  best-fit  disk  model  by  performing  a 
simultaneous  fit  to  the  photometric  and  interferometric 
data.  The  F-star*  s  diameter  and  limb  darkening  are  held 
constant,  aj  and  i  are  free,  whereas  fl,  a;,  e7  P,  and  T  are 
constant. 

5.  Derive  the  best -fit  F-star's  diameter,  limb  darkening 
coefficient,  disk  height,  and  disk  transparency  at  each 
epoch  while  holding  the  remaining  parameters  constant. 
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6.  Finally,  lift  the  constraints  on  our  models  insofar  as 
possible  to  derive  statistical  information  of  the  afore¬ 
mentioned  parameters  via  bootstrapping  individual 
epochs. 


4.  IMAGE  RECONSTRUCTION 

We  performed  image  reconstruction  for  the  figures  presented 
in  this  publication  using  both  SQUEEZE  (Baron  et  al.  2010),  a 
logical  successor  to  the  Markov  Chain  Imager  (MACIM 
Ireland  2006),  and  the  BiSpectrum  Maximum  Entropy  Method 
(BSMEM;  Buscher  1994;  Baron  &  Young  2008).  The  theory 
behind  image  reconstruction,  namely  the  minimization  of  the 
X2  datum  plus  a  regularization  function,  is  common  to  these 
packages;  however,  the  programs  use  different  approaches 
when  solving  the  minimization  problem.  SQUEEZE/MACIM 
perform  global  stochastic  minimization  by  simulated  annealing 
whereas  BSMEM  uses  a  local  gradient-based  approach. 
Despite  the  differences  in  implementation,  the  images  produced 
by  these  packages  are  in  remarkable  agreement.  This  is  proof 
that  given  the  same  data,  the  software  frequently  converge  to 
the  same  solution.  Because  of  this,  we  present  the  SQUEEZE 
total  variation  regularization  images  here  and  other  SQUEEZE 
and  BSMEM  images  as  supplementary  material. 

To  access  the  presence  of  artifacts  in  our  image  reconstruc¬ 
tion  method,  we  adopted  a  pragmatic  approach  in  which  (1)  we 
find  the  best-fit  model  from  our  minimization  process,  (2) 
sample  the  model  using  the  same  UV  coverage  as  the  original 
to  create  synthetic  data,  (3)  redistribute  the  nominal  values  in 
the  synthetic  data  using  the  uncertainties  from  the  real  data,  and 
(4)  reconstruct  the  simulated  model  image  using  the  same 
methods  as  the  real  data.  Steps  two  and  three  were  performed 
using  oifits-simP  The  three  images  (real,  model,  and  simulated 
model)  are  then  compared  qualitatively.  Features  present  in  the 
model  that  are  also  present  in  the  real  and  simulated  model 
images  are  likely  true.  Conversely,  features  seen  in  the  real  or 
simulated  model  image  that  are  not  present  in  the  model  are 
likely  artifacts  of  the  reconstruction  process.  Finally,  features 
present  in  the  real  image  but  not  explained  otherwise  may  be 
real,  but  warrant  further  investigation. 

5.  RESULTS 

5.7.  Model  Selection  for  the  F -star’s  Limb  Darkening  Law 

We  have  used  the  four  post-eclipse  CHARA-MIRC 
observations  to  compare  seven  different  limb  darkening  laws 
for  the  F-star.  As  seen  in  Table  3,  the  uniform  disk  model  is 
always  a  poor  fit  to  the  data  compared  to  the  limb  darkening 
models.  In  all  but  one  epoch,  the  quadratic  limb  darkening  law 
was  found  to  be  a  better  fit  to  the  interferometric  data  than  other 
models.  On  2011  October  10,  the  two-parameter  law  described 
in  Fields  et  al.  (2003)  was  a  better  fit. 

It  would  appear  that  the  quadratic  limb  darkening  law  is 
most  appropriate;  however,  an  inspection  of  the  fits  reveals  that 
all  models  predict  essentially  the  same  visibility  function. 
Interestingly,  all  models  also  predict  significantly  lower  flux  at 
the  limb  than  is  implied  by  either  plane  parallel  or  spherical 
stellar  atmosphere  codes.  We  verified  that  the  SIMTOI  results 
matched  LitPro’s  analytical  models,  therefore,  it  is  unlikely  that 
this  effect  is  fictitious.  We  have  noticed  small,  few  degree,  non- 

13  Oifits  sim:  https://github.com/bkloppenborg/oifits  sim 


zero  closure  phases  and  two  epochs  where  the  location  of  the 
first  visibility  null  differs  between  baselines.  Thus  it  is  possible 
that  the  F-star  may  be  slightly  oblate  or  have  surface  features. 
We  will  explore  these  possibilities  in  greater  detail  in  a  future 
publication.  Because  it  is  simple  and  can  reproduce  the  data, 
we  pragmatically  adopted  the  power-law  limb  darkening  law  to 
represent  the  F-star  in  this  work.  The  mean  out-of-eclipse 
diameter  and  limb  darkening  coefficient  are  2.22  db  0.09  mas 
and  0.50  db  0.26,  respectively.  This  range  of  values  compares 
favorably  with  published  diameters  from  the  NPOI  (#UDD  2.18 
=b  0.05  mas  at  0.5-0.85  /xm,  Nordgren  et  al.  2001),  the  Mark  III 
(limb  darkened  diameter,  #LDD,  1.888-2.136  mas  at 
0.4-0.8  /xm,  Mozurkewich  et  al.  2003),  and  PTI  (#UDD  2.27 
zb  0.11  mas  at  2.2  /xm,  Stencel  et  al.  2008)  interferometers. 

5.2.  Initial  Estimation  of  the  Orbital  Parameters 

Using  the  equations  in  Section  3.4,  we  have  estimated  the 
orbital  semimajor  axis.  Assuming  T  =  9896  zb  1.6  days 
(Stefanik  et  al.  2010),  the  H-band  ingress  time  145  zb  15  days 
(Hopkins  2012),  and  a  circular  orbit,  we  estimate 
olj  =  cq  +  ol2  ~  24  zb  2.6  mas  where  and  a2  are  the 
semimajor  axes  of  the  F-star  and  disk  with  respect  to  the 
system’s  center  of  mass.  In  the  elliptical  case,  the  spectroscopic 
elements  from  Stefanik  et  al.  (2010)  and  Chadima  et  al.  (2010) 
orbital  elements  yield  nearly  identical  results  of  aT  ~  31  zb  3.7 
and  33  zb  4.5  mas,  respectively.  The  dominant  source  of 
uncertainty  in  these  values  comes  from  the  ^10%  errors  in  u 
and  fingress  •  Recognizing  estimates  for  the  F-star’ s  contribution 
to  aT  are  some  13  24  mas  (Strand  1959;  van  de  Kamp  1978; 
Heintz  &  Cantor  1994),  we  establish  bounds  of  13  <  aT  <  38 
mas  for  our  minimization.  A  non-equatorial  intersection,  as 
seen  in  our  previous  work  (Kloppenborg  et  al.  2010),  will 
strictly  decrease  the  upper  bound  on  aT. 

Next,  we  estimated  U  by  fitting  the  in-eclipse  CHARA- 
MIRC  data  with  a  model  consisting  of  a  circular  F-star  with 
power-law  limb  darkening  and  cylindrical  disk.  We  set 
aT  =  31  mas  and  r& Sk  =  10  mas  while  permitting  the  stellar 
and  orbital  inclination  parameters  to  remain  free.  An  average  of 
U  =  296  zb  3  deg  was  obtained  from  the  six  totality  epochs. 

5.3.  Disk  Model  Selection  via  Multi-epoch  Minimization 

From  the  aforementioned  disk  models,  we  created  eight 
variants.  The  geometric  models  are  (1)  a  hard-edged  cylinder, 
and  three  variants  with  power-law  transparency  in  (2)  both 
radius  and  height,  (3)  height  only,  and  (4)  radius  only.  The 
models  from  astrophysical  density  distributions  are  (5)  a 
Pascucci  et  al.  disk  and  (6)  an  Andrews  et  al.  disk.  Although 
not  supported  by  the  interferometric  images,  we  decided  to  test 
for  the  central  clearing  hypothesis  (c.f.  Ferluga  1990,  and 
references  therein)  and  created  (7)  a  Pascucci  et  al.  disk  with  a 
variable  inner  radius.  Finally  we  test  whether  the  Kemp  et  al. 
(1986)  polarization  model  by  creating  (8)  a  Pascucci  et  al.  disk 
that  may  be  tilted  out  of  the  orbital  plane. 

Using  the  MultiNest  minimizer,  we  derived  log  Z  estimates 
by  simultaneously  fitting  each  model  to  a  subset  of  the  H-band 
photometric  data  (consisting  of  200  observations  spaced  at 
approximately  equal  intervals  throughout  the  eclipse)  and  six 
interferometric  epochs  from  CHARA-MIRC  (2009-11,  2009- 
12,  2010-08,  2010-11,  2011-01,  and  2011-09-18).  The  best-fit 
parameters,  posterior  odds  ratio  (A  log  R ),  and  average  x2 
values  for  each  model  are  shown  in  Table  5.  The  A  log  R 
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Table  3 

Posterior  Odds  Ratios  (A  log  R )  Relative  to  the  Uniform  Disk  Model 


Posterior  Odds  Ratio  (A  log  R,  See  Section  3.1) 


Epoch 

Claret  (2000) 

Fields  et  al.  (2003) 

Logarithmic 

Power  Lawa 

Quadratic 

Square  Root 

2011  09  18 

9504 

9518 

9515 

9506 

9545 

9506 

2011  09  24 

1954 

1949 

1959 

1956 

1971 

1954 

2011  10  10 

2568 

2615 

2560 

2561 

2560 

2558 

2011  11  03 

6151 

5900 

6165 

6166 

6227 

6134 

Note. 

a  Although  there  is  slight  evidence  in  favor  of  the  quadratic  limb  darkening  law,  visual  inspection  of  the  visibility  data  shows  no  significant  difference  between  it  and 
the  power  law  limb  darkening  that  we  pragmatically  adopted  in  this  work. 


values  indicate  that  model  8  (the  tilted  Pascucci  disk)  provides 
the  best  simultaneous  fit  to  the  data  and  was  therefore  adopted 
for  the  remainder  of  our  work. 

We  have  rendered  the  best-fit  version  of  each  disk  model  in 
Figure  2.  With  the  exception  of  the  cylindrical  disk,  all  disk 
models  have  a  characteristic  size  of  ^6  mas  radius  and 
~0.75  mas  height  before  becoming  optically  thin.  Plots  of  the 
corresponding  H-band  photometry  in  Figure  3  show  that 
although  the  disk  models  appear  physically  different,  they  all 
do  a  reasonable  job  reproducing  the  global  properties  of  the 
light  curve.  This  clearly  demonstrates  that  reproducing  the 
photometry  is  a  necessary,  but  not  sufficient  condition  to 
prove  the  validity  of  any  particular  disk  model.  We  note  that 
none  of  these  symmetric  disk  models  are  capable  of 
reproducing  all  of  the  photometric  features  seen  during  the 
eclipse. 

Using  model  7  (a  Pascucci  disk  with  a  variable  inner  radius), 
we  test  the  notion  that  the  disk’s  central  clearing  is  responsible 
for  the  alleged  mid-eclipse  brightening.  We  find  that  the  disk 
remains  edge-on,  but  has  an  inner  radius  of  some  3.8  mas.  This 
implies  that  the  inner  60%  of  the  disk  could  be  devoid  of  any 
opaque  material.  Despite  this  fact,  the  impact  on  the  light  curve 
is  minimal.  Given  the  geometry  of  the  eclipse,  light  penetrating 
the  central  clearing  cannot  be  responsible  for  any  mid-eclipse 
brightening.  We  will  discuss  this  result  in  greater  detail  in 
Section  6. 

5.4.  Bootstrapping  and  Aggregate  Statistics 

The  photometry  predicted  by  the  symmetric  disk  models 
creates  an  interesting  corollary:  if  one  assumes  the  disk  is 
symmetric,  one  may  immediately  conclude,  to  the  contrary,  that 
the  disk  must  be  asymmetric  because  the  residuals  between  the 
observed  and  predicted  photometry  far  exceed  the  F-star’s 
AH  rsj  0.05  mag  variations  seen  outside  of  eclipse.  We  tested 
this  asymmetric  conjecture  by  performing  several  single-epoch 
MultiNest  minimizations  to  the  in-eclipse  MIRC  data.  We 
created  three  additional  disk  models  which  (a)  forced  the  disk 
to  reside  in  the  orbital  plane,  (b)  tilted  the  disk  out  of  the  orbital 
plane  at  a  fixed  angle,  and  (c)  permitted  both  the  position  and 
inclination  of  the  disk  with  respect  to  the  orbital  plane  to  vary. 
We  find  that  the  disk  is  tilted  out  of  the  orbital  plane  by  less 
than  4°  (1.33  d=  0.67  deg  with  rejection  of  one  outlier)  and  has 
significant  variations  in  structure  (see  Table  9). 

To  derive  statistical  uncertainties  and  simulated  photometry 
we  used  the  best-fit  values  from  the  aforementioned  minimiza¬ 
tions  as  starting  points  and  bootstrapped  each  interferometric 
epoch  10,000  times.  During  each  bootstrap  we  dynamically 
recalibrated  the  data  using  the  uncertainty  distribution  of  the 
calibrator.  Then  we  created  a  new  realization  of  the  data  using 


Table  4 

A  Summary  of  all  Parameters  Used  in  the  Modeling  Process 


Parameter 

Range 

Units 

Description 

F  star  parameters 

$UDD 

1  3 

mas 

Uniform  disk  diameter 

#LDD 

1  3 

mas 

Limb  darkened  disk  diameter 

GLDD 

0  1 

Power  law  limb  darkening 

coefficient 

Geometric  disk  models 

^in 

0 

mas 

Inner  radius 

fi)Ut 

1  30 

mas 

Outer  radius 

^disk 

0  5 

mas 

Height 

a 

5  20 

Radial  exponent 

0 

0.1  5 

Height  exponent 

Astrophysical  density  disk  models 

KoPo 

3000  10,000 

Characteristic  opacity 

a 

5  20 

Radial  exponent  (Pascucci 

model) 

P 

0.1  5 

Height  exponent  (Pascucci 

model) 

7 

0.001  20 

Radial  exponent  (Andrews 

model) 

K 

0.001  20 

mas 

Disk  scale  height 

rc 

1  4 

mas 

Disk  scale  radius 

ktisk 

±10 

deg 

Disk  inclination 

f^disk 

±10 

deg 

Disk  position  angle 

Orbital  parametersa 

e 

0.227  ±  0.01 

Eccentricy 

i 

70  110 

deg 

Inclination 

id 

39.2  ±  3.4 

deg 

Longitude  of  periastron 

P 

9,896  ±  1.6 

day 

Period 

n 

90  145  and  270  325 

deg 

Position  Angle 

dj,  cq,  d2 

13  <  olj  <  38 

mas 

Semi  major  axis  (total, 

F  star,  disk) 

T 

2,434,723  ±  80 

days 

Time  of  periastron 

Note. 

a  e,  l j,  P,  and  T  from  Stefanik  et  al.  (2010). 


the  measured  uncertainties,  taking  into  account  any  known 
correlations  (i.e.,  as  seen  in  spectrally  dispersed  visibilities)  in 
the  data  when  required.  The  results  of  this  effort  are  shown  in 
Table  6  with  a  subset  of  the  results  plotted  in  Figure  4. 

By  inspection  of  Table  6,  one  can  see  that  the  angular 
diameter  and  limb  darkening  profile  of  the  F-star  are  largely 
consistent  within  lcr.  This  implies  that  there  are  no  egregious 
systematic  calibration  errors  between  our  data  sets.  Further¬ 
more,  these  results  show  that  there  has  been  no  secular  change 
to  the  F-star’s  diameter  over  the  last  14  years.  However,  our 
observations  are  not  sufficiently  precise  to  definitively  exclude 
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Table  5 

Bayes  Factors  Relative  to  the  Cylinder  Model  and  Average  Reduced  x2  for  the  Six  Disk  Models  Described  in  the  Section  3.3 


Name 

Modela 

Orbit 

Disk 

Fit  Information 

OLj 

(mas) 

i 

(deg) 

nn 

(mas) 

r out 

(mas) 

h 

(mas) 

a 

0 

rc 

(mas) 

hc 

(mas) 

Kp  b 

^disk 

(deg) 

fdisk 

(deg) 

A  log  R 

Xr  (H) 

X2r(V2) 

X2,(T3  a) 

xJ(A<,) 

Cylinder 

1 

21.7 

88.4 

0.0 

4.85 

0.54 

0 

260 

10 

12 

48 

Ringed  Disk 

2 

27.0 

88.7 

0.0 

7.20 

0.71 

0.13 

1.64 

63940 

69 

4.1 

4.3 

31 

RingedDisk  (only  (5) 

3 

27.5 

88.8 

0.0 

6.39 

0.79 

0.02 

59811 

93 

4.4 

4.1 

32 

RingedDisk  (only  a) 

4 

27.5 

88.8 

0.0 

7.40 

0.70 

0.11 

63792 

67 

4.1 

4.4 

31 

Pascucci  Disk 

5 

33.0 

89.0 

0.0 

11.11 

1.75 

2.39 

0.029 

6496 

67058 

50 

3.8 

4.1 

30 

Andrews  Disk 

6 

33.0 

89.0 

0.0 

11.14 

1.74 

2.48 

0.032 

5161 

67110 

49 

3.8 

4.1 

31 

Pascucci  Disk  w/  clearing 

7 

32.6 

88.9 

3.8 

10.94 

0.75 

2.32 

0.079 

6667 

67254 

51 

3.8 

4.0 

31 

Tilted  Pascucci  Diskc 

8 

31.2 

88.9 

0.0 

13.33 

3.69 

2.77 

0.007 

6287 

0.02 

2.98 

68806 

50 

3.7 

3.7 

28 

Notes. 

a  See  Section  5.3  for  model  descriptions. 

b  The  product  up  was  treated  as  a  single  quantity  in  our  minimization  process,  see  Section  5.3  for  further  details. 
c  Model  8  obtains  the  highest  Bayes  factor  and  is  therefore  adopted  in  this  work. 
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the  0.6%  year  1  contraction  rate  suggested  by  Saito  & 
Kitamura  (1986). 

In  the  top  panel  of  Figure  4  we  plot  the  H-band  photometry 
as  observed,  simulated  from  the  symmetric  model,  and 
predicted  from  the  bootstrap  process  described  above.  Not  all 
of  the  photometric  values  are  in  perfect  agreement,  but  this  is 
expected  as  the  photometry  was  not  used  as  a  constraint  in  the 
bootstrapping  process.  We  note  that  a  small  change  in  disk 
structure  can  have  a  substantial  impact  on  the  observed 
photometry  (e.g.,  a  0.1  mas  difference  in  thickness,  five  times 
smaller  than  our  resolution  limit,  results  in  AH  ^0.1  mag), 
thus  the  predicted  photometry  is  in  reasonable  agreement. 

In  the  bottom  three  panels  of  Figure  4,  we  plot  the  angular 
diameter,  limb  darkening  coefficient,  disk  scale  height,  hc ,  and 
height  power-law  exponent,  /3,  as  a  function  of  time.  The  1<t 
estimates  for  hc  and  /3  from  model  8  do  not  overlap  well  with 
the  single-epoch  bootstraps.  This  is  likely  due  to  the  ^1.6° 
difference  in  disk  inclinations  between  the  models.  The 
changes  in  hc  and  /3  (greater  hc,  smaller  (3 )  indicate  that  the 
disk  is  more  spatially  extended  in  height  before  mid-eclipse 
than  it  is  after  mid-eclipse.  Combined  with  the  asymmetric 
evolution  of  several  neutral  absorption  lines  during  eclipse  (c.f. 
Lambert  &  Sawyer  1986;  Leadbeater  et  al.  2012),  the  evidence 
suggests  that  the  disk  is  not  purely  symmetric.  We  suspect  that 
these  features  could  be  explained  by  asymmetric  heating  (c.f. 
Takeuti  2011)  and  sublimation  of  the  disk  on  the  side  facing  the 
F-star.  Verifying  this  claim  by  hydrodynamical  radiative 
transfer  simulations  is  beyond  the  scope  of  this  work. 

Because  of  the  limited  UV  coverage,  the  egress  CHARA- 
CLIMB  data  was  fit  in  conjunction  with  an  interpolated 
photometric  point.  The  result  reveals  that  the  F-star  is  similar  in 
size,  but  the  disk  is  more  extended  in  the  vertical  direction  than 
the  CHARA-MIRC  observations  three  months  earlier  would 
imply.  This  interpretation  is  supported  by  the  appearance  of  the 
model-independent  images  (discussed  below);  however,  we 
caution  the  reader  that  this  may  simply  be  an  artifact  of  the 
limited  UV  coverage  of  this  data  set. 

In  Table  7,  we  show  the  aggregate  statistics  derived  from 
single-epoch  bootstrapping  for  each  beam  combiner.  We 
believe  these  values  represent  the  general  characteristics  of 
the  system.  In  Table  8,  we  summarize  these  quantities  for  a 
variety  of  distance  estimates  for  the  system.  Due  to  the  large 
scatter  of  possible  distances  (0.5  4kpc),  we  use  only  the 
nominal  values  and  do  not  propagate  any  uncertainties  from  the 
distance  measurements.  We  caution  the  reader  that  the 
aggregation  of  data  to  create  Table  8  was  performed  without 
regard  to  either  the  wavelength  of  observation  or  any 
asymmetries  we  advocate  exist.  Hence,  the  outliers  have  biased 
and  skewed  the  resulting  value.  We  provide  this  last  table  to 
assist  with  the  creation  of  a  full  radiative  transfer  model  of  the 
system  including  dust  physics  rather  than  provide  a  definitive 
measurement  of  the  properties  of  the  system. 

5.5.  Reconstructed  Images  and  Artifacts 

In  Figure  5,  we  present  the  best-fit  model  and  model- 
independent  SQUEEZE  reconstructions  using  the  total  varia¬ 
tion  regularizer.  By  following  the  qualitative  comparison 
method  discussed  in  Section  4,  we  may  distinguish  whether 
the  features  seen  in  these  images  are  real  signals  or  artifacts  of 
the  reconstruction  process.  A  detailed  discussion  of  this 
process,  as  well  as  BSMEM  and  other  SQUEEZE  reconstruc¬ 
tions,  can  be  found  in  the  Appendix  and  extended  Figures 


A46  A52.  Summarizing  this  account,  the  true  features  are  as 
follows:  (1)  the  dark  lane  in  the  F-star’ s  southern  hemisphere, 
interpreted  to  be  the  disk;  (2)  flux  that  appears  on  the  far  south 
west  (or  south  east)  edge  of  the  F-star  during  ingress  (egress); 
and  (3)  the  presence  or  absence  of  the  F-star’ s  southern  pole. 
Although  the  post-eclipse  interferometric  observations  do 
feature  small  non-zero  closure  phases  and  photometry  of  the 
e  Aurigae  system  does  show  an  intrinsic  AV  ^  0.1  mag 
photometric  variation  outside  of  eclipse,  we  presently  regard 
any  flux  variations  on  the  surface  F-star  as  artifacts. 

6.  CONCLUSIONS  AND  DISCUSSION 

We  have  analyzed  106  nights  of  interferometric  observations 
provided  by  four  beam  combiners  at  three  different  interfero¬ 
metric  facilities  to  derive  properties  of  the  F-star,  determine 
previously  unknown  orbital  elements,  and  access  the  global 
structures  of  the  optically  thick  portion  of  the  eclipsing  disk. 
We  have  reconstructed  a  series  of  model-independent  images 
using  SQUEEZE  and  BSMEM  programs  with  a  variety  of 
regularization  functions.  The  images  show  that  the  F-star 
appears  circular  (round)  and  free  of  egregious  asymmetries 
prior  to  and  after  the  2009  2011  eclipse.  During  the  eclipse, 
most  of  the  southern  hemisphere  on  the  F-star  is  obscured.  The 
appearance  and  persistence  of  this  feature  lead  us  to  confirm  the 
Huang  Lissauer  disk  hypothesis  for  the  eclipse.  During  most  of 
the  epochs,  the  southern  pole  of  the  F-star  is  visible,  thereby 
providing  an  opportunity  to  measure  the  thickness  of  the  disk  at 
sub-milliarcsecond  resolution. 

Our  interferometric  modeling  efforts  were  complex:  under 
a  Bayesian  framework,  we  tested  the  observations  of  the 
F-star  against  seven  different  analytic  limb  darkening 
prescriptions  and  differentiated  between  eight  proposed  disk 
models.  The  Bayes  factors  listed  in  Table  5  are  all 
exceptionally  large.  After  conducting  a  comprehensive  over¬ 
view  of  SIMTOV  s  rendering  pipeline,  liboV  s  unit  test 
framework,  and  our  use  of  MultiNest,  we  found  no  mistakes 
in  our  implementation.  Hence  these  values  are  either  true,  or 
there  is  some  unforeseen  systematic  error  for  which  or 
modeling  process  did  not  account.  Nevertheless,  we  are 
confident  that  our  best-fit  model  is  indeed  the  most  probable 
as  it  achieves  the  lowest  reduced  chi-squared  estimate.  The 
pre-  and  post-eclipse  observations  indicate  that  the  F-star  has 
remained  at  a  more-or-less  constant  diameter  for  the  last  14 
years;  however,  our  observations  are  not  sufficiently  precise 
to  definitively  exclude  the  0.6%  year  1  contraction  rate 
suggested  by  Saito  &  Kitamura  (1986).  The  average  power- 
law  limb  darkening  coefficient,  ^0.5,  is  much  higher  than 
predicted  for  an  F-type  supergiant  (^0. 1-0.2).  The  presence 
of  small,  <5°,  non-zero  closure  phases  on  the  longest 
baselines  in  the  post-eclipse  observations,  coupled  with 
variations  in  radius  and  limb  darkening  seen  during  the 
eclipse,  suggest  that  there  may  be  convective  cells  or  some 
other  feature  (e.g.,  spots)  on  the  surface  of  the  F-star.  These 
features,  if  they  exist,  have  only  a  minimal  impact  on  our 
results. 

We  have  simultaneously  fit  the  disk  models  to  a  subset  of  the 
interferometric  and  photometric  data.  We  find  that  the  data  can 
be  adequately  fit  by  a  variety  of  models;  however,  the  most 
consistent  model  was  that  of  tilted  disk  derived  from  an 
astrophysical  density  distribution.  The  opaque  region  of  the 
disk  is  seen  nearly  edge  on  and  is  remarkably  uniform.  These 
conclusions  are  supported  by  spectroscopic  work  by  Griffin  & 
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Figure  2.  (Top)  The  eight  best  fit  symmetric  disk  models  resulting  from  a 
simultaneous  minimization  to  both  the  photometric  and  interferometric  data. 
The  models  are  as  described  in  Section  3.3.  All  models  have  been  rotated  in 
position  angle,  but  otherwise  appear  as  they  would  when  occulting  the  F  star's 
photosphere.  Model  #8  has  the  greatest  evidence  value  and  was  adopted  for 
the  remainder  of  this  work.  (Bottom)  The  disk's  silhouette  created  by  aligning 
the  individual  SQUEEZE  reconstructions  with  our  orbital  solution.  See 
Appendix  B  for  further  details. 


Date  (JD  -  2450000.0) 


Figure  3.  Observed  and  simulated  H  band  photometry  from  the  eight  best  fit 
symmetric  disk  models  as  a  function  of  time.  The  four  rows  of  tick  marks  at  the 
bottom  of  the  plot  indicate  the  time  of  observation  by  the  corresponding  beam 
combiner.  Tick  marks  appear  thicker  when  successive  nights  of  observations 
occurred.  None  of  the  symmetric  disk  models  entirely  reproduce  all  features  of 
the  light  curve,  implying  that  the  disk  is  not  symmetric.  The  F  star's 
AH  ™  0.05  mag  out  of  eclipse  variations  are  clearly  seen  al  all  phases  of  the 
eclipse. 


Stencel  (2013),  who  have  shown  the  precise  repetition  of  disk- 
related  spectral  features  have  not  changed  appreciably  over  the 
last  century.  Therefore,  it  is  unlikely  that  the  disk  is 
significantly  twisted  or  warped  (Kumar  1987).  We  do, 
however,  find  evidence  that  the  disk  may  be  slightly  tilted 
out  of  the  orbital  plane.  If  this  is  true,  the  difference  between 


our  model  and  the  Kemp  et  al.  (1986)  polarization  result  could 
be  attributed  to  precession. 

The  thickness  and  inclination  of  the  disk  exclude  the 
possibility  that  a  the  purported  mid-eclipse  brightening  is 
caused  by  light  penetrating  a  central  clearing  in  the  disk  (cf. 
Wilson  1971)  or  the  notion  that  a  series  of  semi-transparent 
rings  are  responsible  for  the  photometric  variations  seen  during 
totality  (cf.  Ferluga  1990,  1989).  We  suggest  that  these  light 
curve  features,  if  true,  have  other  physical  causes.  For  example, 
a  mid-eclipse  brightening  could  be  due  to  scattering  above  and 
below  the  plane  of  the  disk  (e.g.,  Budaj  201 1;  Muthumariappan 
&  Parthasarathy  2012),  perhaps  in  the  same  region  responsible 
for  the  increase  in  He  10830  A  absorption  (Stencel  et  al.  201 1). 
Likewise,  the  manifestation  of  AV  ^0.1  mag  variations  during 
totality  are  probably  orbitally  excited  non-radial  pulsation  of 
the  F-star  (Kloppenborg  et  al.  2012),  rather  than  substructure  in 
the  disk. 

Finally  we  predict  that  the  secondary  eclipse  will  occur 
between  ^JD  2,461,030—2,461,860  (2025  December  20  2028 
March  29).  We  encourage  a  comprehensive  photometric 
campaign  during  this  time  focusing  on  NIR,  mid-IR,  and  far- 
IR  observations  to  confirm  this  prediction. 

The  CHARA  Array,  operated  by  Georgia  State  University, 
was  built  with  funding  provided  by  the  National  Science 
Foundation  grant  AST -0606958,  Georgia  State  University, 
the  W.  M.  Keck  Foundation,  and  the  David  and  Lucile 
Packard  Foundation.  This  research  is  supported  by  the 
National  Science  Foundation  as  well  as  by  funding  from 
the  office  of  the  Dean  of  the  College  of  Arts  and  Science  at 
Georgia  State  University.  MIRC  was  supported  by  the 
National  Science  Foundation.  The  Navy  Precision  Optical 
Interferometer  is  a  joint  project  of  the  Naval  Research 
Laboratory  and  the  US  Naval  Observatory,  in  cooperation 
with  Lowell  Observatory,  and  is  funded  by  the  Office  of 
Naval  Research  and  the  Oceanographer  of  the  Navy. 
Participants  from  the  University  of  Denver  are  grateful  for 
the  bequest  of  William  Hershel  Womble  in  support  of 
astronomy  at  the  University  of  Denver.  They  acknowledge 
support  from  National  Science  Foundation  through  1SE  grant 
DRL-0840188  to  the  American  Association  of  Variable  Star 
Observers  and  AST  grant  10-16678  to  the  University  of 
Denver. 

Facilities:  CHARA,  NPOl,  PQ:PT1. 

APPENDIX  A 

IMAGE  RECONSTRUCTION  AND  ARTIFACT 
DISCUSSION 

In  the  following  section,  we  elaborate  on  the  image 
reconstruction  and  artifact  detection  process  for  all  epochs. 
We  present  the  best-fit  S1MTOI  model  and  reconstructed  with 
SQUEEZE  and  BSMEM.  All  images  are  rendered  with  North 
up  and  east  to  the  left.  The  0.5  mas  H-band  or  0.7  mas  K-band 
resolution  limit  of  the  data  is  indicated  by  the  circle  in  the  lower 
left-hand  comer  of  the  model  image. 

The  SQUEEZE  reconstructions  were  conducted  using  the 
Laplacian  (abbreviated  LA),  Total  Variation  (abbreviated 
TV),  and  L0-norm  (abbreviated  L0)  regularizes.  Images 
were  reconstructed  in  a  grid  with  regularizer  weights  ranging 
from  0.1  100,000  in  semi- logarithmic  steps.  The  images  with 
the  lowest  reduced  chi-squared  were  selected  for  presenta¬ 
tion.  Most  of  the  epochs  were  reconstructed  from  a  Dirac 
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Table  6 

Bootstrapped  Nominal  Values  and  Uncertainties  for  Model  8  With  a  Fixed  Tilt  Evaluated  on  a  Per  epoch  Basis  Subject  to  Only  the  Interferometric  Data 


F-star 

Disk 

Statistical  Information 

Data  set 

V(V2)  N(T3)  N(UV ) 

Effective  JD 

#UDD 

(mas) 

#ldd  oldd  kP 

(mas) 

(  x  IQ"2)  (  x  lO-2) 

rc 

(mas) 

K 

(mas) 

a 

X2r 

X2r(V2)  X2(Tm)  x\(T^)  Notes 

1997-10- 

14 

16 

2450744.0200 

2.43 

0.38 

0.38 

22-PTI 

± 

0.29 

1998-11- 

66 

69 

2450761.9500 

2.13 

0.27 

0.27 

Distribution  not  well 

07-PTI 

± 

constrained 

0.28 

1997-11- 

20 

21 

2451124.9500 

2.78 

2.21 

2.21 

09-PTI 

± 

0.15 

1998-11- 

10 

11 

2451142.9400 

1.93 

0.78 

0.78 

Distribution  well  con¬ 

25-PTI 

m 

strained.  HD  30823 

0.44 

sole  calibrator  this 
night.  Perhaps  cali¬ 
brator  diameter 

over-estimated? 

1998-11- 

5 

5 

2451143.9300 

2.06 

0.53 

0.53 

Distriubtion  poorly 

26-PTI 

± 

constrained.  Nom¬ 

0.53 

inal  value  matches 

best-fit  MultiNest 

estimate. 

2006- 

540 

135 

544 

2453791.7434 

2.09 

2.35 

2.65 

2.00 

1.49 

02-NPOI 

± 

0.06 

2007- 

660 

330 

966 

2454173.6316 

9  OQ+0.07  rt  /19+O.23 

Z.25_o  o2  t).4Z_0  05 

2.28 

0.90 

1.67 

5.67 

03-NPOI 

2007-10- 

100 

104 

2454392.9900 

2.13 

0.59 

0.59 

19-PTI 

± 

0.13 

2007-10- 

81 

84 

2454393.9700 

2.08 

0.18 

0.18 

20-PTI 

± 

0.14 

2007-10- 

40 

43 

2454394.9800 

2.16 

1.35 

1.35 

21-PTI 

i 

0.24 

2007-11- 

20 

22 

2454431.8500 

2.55 

2.16 

2.16 

Distribution  not  con¬ 

27-PTI 

± 

strained,  highly 

0.32 

skewed  toward 
higher  values. 

2007-12- 

5 

5 

2454457.7800 

2.09 

0.68 

0.68 

23-PTI 

± 

0.43 

2007-12- 

35 

39 

2454458.7500 

2.45 

1.75 

1.75 

24-PTI 

± 

0.26 

2008-02- 

10 

10 

2454513.6400 

1.75 

3.22 

3.22 

Distribution  is  well 

17-PTI 

± 

constrained.  Bad 

0.33 

calibration? 

2008-02- 

25 

27 

2454514.6500 

2.17 

1.77 

1.77 

18-PTI  ± 

0.34 


H 

X 


Z 

H 

C/3 

W 

E 

w 


Table  6 

(Continued) 


Data  set 

A(V2) 

N(T3) 

N(UV) 

Effective  JD 

2008-09- 

CHARA- 

MIRC 

23 

8 

50 

2454729.0153 

2008-10- 

17-PTI 

65 

69 

2454757.0000 

2008-10- 

26-PTI 

80 

84 

2454765.9800 

2008-11- 

CHARA- 

MIRC 

138 

76 

268 

2454778.5967 

2008-11- 

08-PTI 

10 

11 

2454778.8700 

2008-11- 

09-PTI 

20 

21 

2454779.8700 

2008-11- 

16-PTI 

10 

10 

2454786.8700 

2008-11- 

22-PTI 

50 

52 

2454792.8200 

2008-12- 

CHARA- 

MIRC 

38 

8 

65 

2454810.8249 

2009- 

03-NPOI 

840 

420 

1225 

2454898.6326 

2009-11- 

CHARA- 

MIRC 

1091 

672 

2575 

2455138.9326 

2009-12- 

CHARA- 

MIRC 

730 

392 

1662 

2455169.0375 

2009- 

12-NPOI 

290 

29 

293 

2455185.6210 

2010- 

01  AB- 

NPOI 

3324 

1376 

4444 

2455205.5082 

2010- 

02-NPOI 

810 

265 

814 

2455243.1816 

2010-02- 

CHARA- 

MIRC 

96 

64 

236 

2455245.7444 

#UDD 

(mas) 


2.37 

± 

0.20 

2.01 

i 

0.15 


2.17 

i 

0.40 

2.34 

± 

0.49 


2.30 

db 

0.23 


F-star  Disk 


#ldd  oLDD  up  rc  hc 

(mas)  (mas)  (mas) 

_ (  x  10~2)  (  x  10~2) _ 

2.2S±8;81  0.69±8;!8 


2.22±8.06  0-39±8:iI 


2-36±8;8I  0.80±g;J? 


9  1  6+0-07  0  97+°-25 

z-lo-0.02  '“'••2  '—0.04 


9  9Q+0.02 
z-zy-0.03 

0.62^ 

6329^2038 

i  7Q+0.02 
A-/y-0.02 

9  C7+0.40 
z-9 '  — 0.40 

9-19±g;8t 

9  9Q+0-04 

+.+o_o.o4 

l-00+g;B 

6842^2009 

9  79+0-07 

3-lz-0.01 

12. 18^i  oo 

19.94lg;S 

9  ns+0-21 
Z.U9_o.04 

0.72^ 

7007^}§47 

9  99+0-22 

9-z9_o.03 

6.70±l.fg 

16.72+88 

9  94+0-03 

o.o2 

0.76^0.07 

621li!Io° 

9  Q7+0.08 
z-y/-0.06 

4.io+g|g 

17.84±8;gl 

9  97+0-10 

z-9'-0.10 

0.62 ±g;U 

62191S 

io.7o±;j8 

2.oii§a 

o.2iig:l| 

3075±}ggg 

13. 14+8.i8 

Statistical  Information 

+  x2r  X2r(V2)  X2Ma)  X2r(T34,)  Notes 


1.91 

2.46 

1.94 

0.29 

Errors  limited  by  cali¬ 
brator  uncertainty 

5.71 

5.71 

0.88 

0.88 

3.81 

4.59 

3.24 

2.96 

Errors  limited  by  cali¬ 
brator  uncertainty 

0.08 

0.08 

2.18 

2.18 

Distribution  not  con¬ 
strained  and  highly 
skewed  toward 
higher  values.  Bad 
calibration? 

1.32 

1.32 

3.67 

4.03 

5.42 

0.23 

Errors  limited  by  cali¬ 
brator  uncertainty 

1.85 

1.65 

1.29 

2.80 

i.56±g:;§ 

4.01 

2.64 

0.00 

8.57 

o.io±8;§! 

9.26 

4.08 

3.77 

24.39 

0.92 

2.14 

2.20 

1.89 

1.73 

Visibities  at  short 

baselines  are  much 
higher  than  CHARA 
model  would 
predict. 


2.43^$ 

1.23 

1.61 

1.21 

0.33 

Disk  rc  is  clearly 
bimodal 

0.70+023 

1.42 

1.47 

1.81 

0.87 

0.55tg;!i 

6.88 

1.64 

1.18 

20.45 

Interferometry  +  Pho- 

tometry.  Statistics 
from  MultiNest 
distribution. 
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Data  set 


2010- 

04-NPOI 

2010-08- 

CHARA- 

MIRC 

2010- 09- 
CHARA- 
MIRC 

2010-10- 

CHARA- 

MIRC 

2010-11- 

CHARA- 

MIRC 

2010-12- 

CHARA- 

MIRC 

2011- 01- 
CHARA- 
MIRC 

2011-04- 

CHARA- 

CLIMB 

2011-09-18- 

CHARA- 

MIRC 

2011-09-24- 

CHARA- 

MIRC 

2011-10-10- 

CHARA- 

MIRC 

2011-11-03- 

CHARA- 

MIRC 


Table  6 

(Continued) 


F-star 

Disk 

Statistical  Information 

N(V2) 

N(T3) 

N(UV) 

Effective  JD 

#UDD 

#LDD 

<+LDD 

K,p 

rc 

K 

(? 

x2r(v2) 

x)(TM) 

X)(T^)  Notes 

(mas) 

(mas) 

(  X  10"2) 

(  X  10"2) 

(mas) 

(mas) 

15 

0 

15 

2455289.6026 

9  99+0.28 
z-4-5-0.28 

0.56±§J} 

6500+?ggg 

n.5oll;lg 

0.8811:5! 

0.53 

0.53 

0.00 

0.00 

960 

640 

2164 

2455430.5170 

9  9  9+0.04 
z-J:)-0.04 

0.74ig;?? 

99oii!!|J 

1  97+O.IO 
A-Z/-0.10 

9  7  r +0.07 

4-  '4—0.12 

10.74 

3.97 

3.90 

21.1  A 

1176 

728 

3020 

2455464.4883 

9  97+O.O3 
Z-J  '—0.02 

0.73ig;§S 

59441IH! 

9  9C+O.2O 
z-3-3— 0.10 

9  1  1+0.10 

4-11-o.os 

3.38 

1.96 

1.98 

7.07 

288 

152 

732 

2455496.4319 

2.26±g;gl 

0.62±g;gl 

6629i{997733 

1  4.7+O.8O 

3.56lgil 

1.77 

1.31 

0.77 

3.67 

288 

192 

763 

2455505.4193 

9  i  o+O.Ol 
Z.l5_o.oi 

1  1+0.05 
u- 0.03 

372011H? 

9.6618:18 

o.4o15:J? 

4.64 

2.93 

1.16 

10.71 

191 

112 

475 

2455543.7059 

9  99+0.04 

z-zz— 0.03 

O  9 1+0-10 
H-3 1-0.09 

5548llg| 

0.7218:18 

4.7118.1! 

5.72 

2.29 

1.36 

15.92 

310 

182 

860 

2455580.2465 

9  i  o+0.02 
Z.15_002 

n  9C+O.O4 
u-z:)-0.04 

301311!!! 

4.94lJig 

1.84181! 

5.01 

2.57 

1.36 

12.83 

41 

14 

45 

2455655.0673 

9  99+0.06 

z-42>-0.06 

0.33;gj| 

elooliggg 

9  a  O+0.30 

4  Ao_ 0.30 

13.20ll.l8 

18.55111!  0.64lg  ?g 

3.23 

4.06 

1.04 

2.97 

201 

240 

756 

2455823.0305 

9  9C+0.02 
z,z-3— 0.04 

0.62±g;?§ 

5.06 

3.55 

3.02 

8.37 

120 

160 

394 

2455829.0277 

9  1 7+0.03 
Z1  7 -0.03 

0.36±g$ 

5.40 

3.32 

1.67 

10.69 

400 

480 

1412 

2455844.9422 

9  1  9+0.04 
Z-iZ-0.05 

0.34 

16.97 

5.49 

6.82 

36.70 

831 

1119 

2677 

2455868.8509 

9  9C+O.O3 
z,z  — 0.07 

0.57±g;?f 

12.14 

5.64 

7.21 

21.90 

Note. 

a  The  changes  seen  in  height  power,  j3,  hint  that  there  may  be  some  asymmetric  structure  in  the  disk. 
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Date  (JD  -  2450000.0) 

Figure  4.  Results  of  single  epoch  minimizations  to  the  interferometric  (and  in 
two  cases  photometric)  data.  Top:  the  observed  and  predicted  H  band 
photometry  as  a  function  of  time.  Superimposed  is  the  best  fit  symmetric  disk 
model  (#8)  and  predicted  H  band  magnitudes  from  the  single  epoch  fits. 
Bottom:  the  best  fit  F  star  diameter  (limb  darkened  or  uniform);  limb  darkening 
coefficient;  disk  scale  height;  and  disk  radial  exponent,  ft,  as  a  function  of  time. 
The  horizontal  bands  denote  the  la  limits  on  the  parameters  from  model  eight. 
The  diameter  of  the  F  star  is  largely  consistent,  regardless  of  the  wavelength  of 
observation.  This  indicates  there  is  not  an  egregious  systematic  calibration  error 
between  the  different  interferometric  facilities. 

starting  image.  However,  the  sparse  UV  coverage  in  the 
2008-11,  2010-02,  and  201 1-04  mandated  we  use  the  best-  fit 
model  images  to  initialize  the  flux  distribution.  SQUEEZE 
was  executed  with  50  threads  with  500  realizations  each. 

All  BSMEM  reconstructions  were  performed  using  flat 
priors  and  a  1.4  mas  diameter  Gaussian  for  entropy  estimation. 
We  have  used  the  “full”  elliptical  approximation  for  the 
bi spectra  uncertainties  (see  BSMEM  documentation  for  details) 
which  we  found  dramatically  suppress  reconstruction  artifacts. 

A.L  2008-1 1 

Figure  A46  shows  the  resulting  model  and  image  from  this 
set  of  two  four-tele  scope  MJRC  observations  that  were  taken 
at  nearly  the  same  hour  angle.  The  image  reconstruction  in 
this  epoch  is  quite  poor  when  compared  with  later  epochs.  By 
inspection  of  the  SED  (Hoard  et  al.  2012,  2010)  we  know  that 
the  F-star  is  the  dominant  source  of  flux  in  the  H  band,  hence 
we  interpret  this  to  be  an  image  of  the  F-star.  All  of  the 
images  show  that  the  flux  is  mostly  constrained  within  the 


bounds  of  the  best-fit  model.  The  model-independent  image 
shows  the  star  as  approximately  round  with  some  surface 
features.  These  features  are  mostly  replicated  in  the  model- 
independent  reconstruction  from  the  synthetic  data,  hence  the 
non -circular  structure  and  small  photometric  variations  seen 
on  the  F-star  are  artifacts  of  the  reconstruction  process. 
Therefore  the  F-star  before  the  eclipse  does  not  appear  to 
have  any  egregious  asymmetries  which  might  interfere  with 
late  r  ob  ser vati  on  s . 

A. 2.  2009-11 

In  Figure  A47,  we  show  the  best- fit  model  and  model- 
independent  images  reconstructed  from  this  eclipse  ingress 
phase  epoch.  For  reasons  discussed  above  (and  in  greater  detail 
in  Kloppenborg  et  al.  2010),  we  again  interpret  the  bright 
source  to  be  the  F-star  and  the  dark  region  occuring  in  the 
southern  half  of  the  image  as  the  disk  intruding  into  the  line  of 
sight.  We  are  not  aware  of  any  evidence  which  suggests  the 
F-stars  rotation  is  misaligned  to  the  binary's  orbit,  hence  we 
shall  call  the  un-eclipsed  portion  of  the  F-star  the  “northern 
hemisphere.”  Likewise,  we  will  refer  to  the  (mostly)  eclipsed 
portion  as  the  “southern  hemisphere.”  The  North  pole  of  the 
F-star  would  be  located  at  a  position  angle  of  ^26°. 

By  comparing  the  best -fit  model  from  S1MT01  to  the 
sampled  model  reconstruction,  we  may  qualitatively  access  the 
presence  of  artifacts  in  the  image.  The  dark  spot  in  the  northern 
hemisphere  and  two  bright  spots  in  the  east/ west  near  the  limb 
appear  to  be  artifacts.  The  straight  edges  along  the  perimeter  of 
the  F-star  are  a  common  artifact  caused  by  the  UV  coverage  of 
the  data  set. 

Despite  the  large  number  of  artifacts,  several  real  features 
may  be  discerned.  For  example,  the  southern  pole  is  seen  in  the 
real  image,  model,  and  synthetic  reconstruction;  hence  we  feel 
this  feature  is  real.  Likewise,  the  small  amount  of  flux  seen  on 
the  western  edge  of  the  disk  intrusion  is  also  real. 

AJ.  2009-12 

The  best-fit  model  and  reconstructed  images  of  this  second 
ingress  epcoh  are  shown  in  Figure  A48.  This  data  has  excellent 
UV  coverage  and  appears  similar,  in  many  regards,  to  the 
2009-11  epoch.  Like  previous  observations,  the  straight-edge 
appearance  of  the  F-star  is  an  artifact  of  the  UV  coverage.  It  is 
probable  that  the  spots  seen  in  the  northern  hemisphere  of  the 
F-star  are  also  artifacts. 

The  southern  pole  again  appears  quite  strong  in  the  real 
image,  model,  and  synthetic  reconstruction,  implying  this 
feature  is  likely  real.  For  the  same  reason,  we  regard  the  small 
quantity  of  flux  at  the  west  edge  of  the  disk  intrusion  to  be  a 
real  feature  rather  than  an  artifact  of  the  reconstruction  process. 

A.4.  20104)2 

The  UV  coverage  at  this  epoch  is  extremely  poor,  consisting 
of  two  four-telescope  observation  with  M1RC.  Hence  the 
model-independent  images  shown  in  Figure  A49  are  difficult  to 
interpret  without  information  garnered  from  the  model  and 
H-band  photometry.  Both  the  model  and  images  imply  that  the 
entire  southern  hemisphere  and  a  small  fraction  of  the  northern 
hemisphere  are  covered.  This  conclusion  is  supported,  at  least 
circumstantially,  by  the  H-band  photometry  being  at  its  faintest 
at  this  time.  The  appearance  of  bright  spots  in  the  northern 
hemisphere  is  most  likely  caused  by  limited  UV  coverage  or 
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Table  7 

Aggregate  Statistics  for  all  Interferometric  Data  with  Uncertainties  Determined  from  the  Maximum  of  the  Upper/Lower 
Averaged  Bootstrapped  Uncertainties  or  the  Standard  Deviation  of  the  Nominal  Values 


Quantity 

Units 

NPOI  (V) 

MIRC  (H) 

CLIMB  (K) 

PTI  (K) 

Quantity 

NPOI  (V) 

MIRC  (H) 

CLIMB  (K) 

PTI  (K) 

#UDD 

(mas) 

2.09  ±  0.06 

2.10  ±  0.15 

2.22  ±  0.53 

#LDD 

(mas) 

2.21  ±  0.28 

2.22  ±  0.09 

2.33  ±  0.06 

gldd 

0.47  ±  0.28 

0.50  ±  0.26 

0.33  ±  0.48 

Udisk 

(deg) 

1.30  ±  0.67 

Uisk 

(deg) 

89.49  ±  1.03 

Kpa 

6676  ±  2000 

5667  ±  2188 

6100  ±  2000 

rc 

(mas) 

3.10  ±  0.22 

2.76  ±  1.36 

3.48  ±  0.30 

K 

(mas) 

0.07  ±  0.09 

0.05  ±  0.05 

0.13  ±  0.04 

a 

17.28  ±  0.79 

14.56  ±  7.60 

18.55  ±  1.96 

0 

1.23  ±  1.05 

2.18  ±  1.67 

0.64  ±  0.50 

Notes.  Orbital  values  of  0 

=  297.60  ±  0.06  (deg), 

i  =  88.89  ±  0.03  (deg),  and  aT  =  31.2  ±  0.9  (mas)  were  used  in  these  models. 

a  Not  well  constrained. 

Table  8 

Representative1 

1  Linear  Equivalent  of  our  Results  if  the  System  were  at  Various  Distances  in  Literature 

Nominal  distance  estimates  (pc) 

Quantity 

This  work 

600b,c 

653d  737e 

1000f 

1500s 

Linear  units 

Q 

297  ±  3 

(deg) 

i 

89  ±  1 

(deg) 

olj  =  cq  +  oli 

31  ±  3 

(mas) 

18.72  ±  1.80 

20.37  ±  1.96  22.99  ±  2.21 

31.20  ±  3.00  46.80  ±  4.50 

(AU) 

F  star  Radius 

1.11  ±  0.05 

(mas) 

143.25  ±  5.81 

155.90  ±  6.32  175.96  ±  7.13 

238.75  ±  9.68  358.13  ±  14.52 

(*©) 

F  star  LDD  coeff 

0.50  ±  0.26 

Disk  scale  Height  ( hc ) 

1.038  ±  0.139 

(mas) 

0.03  ±  0.03 

0.03  ±  0.03  0.04  ±  0.04 

0.05  ±  0.05 

0.07  ±  0.07 

(AU) 

Disk  scale  radius  (rc) 

7.416  ±  0.276 

(mas) 

1.66  ±  0.82 

1.80  ±  0.89  2.03  ±  1.00 

2.76  ±  1.36 

4.14  ±  2.04 

(AU) 

Notes.  These  values  average  over  all  interferometric  epochs.  Therefore,  these  estimates  are  biased  and  are  skewed  by  the  outliers  in  Table  6. 
a  We  caution  the  reader  that  this  aggregation  of  data  is  performed  in  a  wavelength  and  model  agnostic  fashion.  Thus,  any  asymmetries  in  the  system,  which  we  argue 
exist,  have  biased  the  values  quoted  here.  These  values  are  supplied  to  ease  the  creation  of  a  radiative  transfer  model  including  dust  physics.  We  do  not  advocate  that 
these  values  be  quoted  elsewhere. 
b  van  de  Kamp  (1978). 
c  Heintz  &  Cantor  (1994). 
d  van  Leeuwen  (2008). 
e  Kloppenborg  (2012). 
f  Strand  (1959). 

8  Guinan  et  al.  (2012). 


the  reconstruction  process,  rather  than  any  real  surface  flux 
variations  on  the  F-star.  The  model  independent  and  synthetic 
images  agree  quite  well  about  the  over-all  appearance  of  the 
F-star  during  this  epoch. 


A.5.  2010-08  /  2010-09  /  2010-10  /  2010-11  /  2010-12  / 
2011-01 

The  qualitative  appearance  of  these  epochs  is  quite  similar 
(see  Figure  A50),  hence  they  will  be  discussed  in  aggregate. 
The  obscuration  by  the  disk  remains  remarkably  consistent 
across  five  months  of  observations.  The  occasional  spot  in  the 
F-stars  northern  hemisphere,  scalloped  edge  of  the  disk  along 
the  F-stars  equator,  and  flux  variations  along  the  F-stars  equator 
are  frequently  seen  in  the  real  data  and  synthetic  reconstruc¬ 
tions,  hence  these  are  likely  artifacts.  The  southern  pole  has  re¬ 
appeared.  It  appears  in  the  real  image,  model,  and  synthetic 
image  therefore  we  regard  this  as  a  true  feature  in  the  image. 

It  is  important  to  note  that  the  mid-eclipse  observation 
(2010-08)  shows  the  disk  as  entirely  opaque.  Hence  mid¬ 
eclipse  brightening  hypotheses,  which  rely  on  a  large  opening 


in  the  disk,  are  unlikely.  Likewise,  although  photometric 
variations  were  seen  during  this  phase  of  the  eclipse,  no 
significant  flux  variations  were  seen  within  the  disk  plane, 
hence  the  photometric  variations  are  not  likely  a  result  of  flux 
penetrating  semi-transparent  gaps  in  the  disk’s  midplane. 

A.6.  2011-04 

In  Figure  A5 1  we  present  the  only  interferometric  images  of 
the  e  Aurigae  during  the  egress  phase  and  the  first  image 
created  with  data  taken  by  the  CLIMB  beam  combiner.  Unlike 
all  other  data  sets,  the  LO-norm  images  also  required  the  use  of 
the  uniform  disk  regularizer  in  SQUEEZE.  These  model- 
independent  K-band  images  show  that  a  portion  of  the  F-star’ s 
south  east  edge  is  no  longer  obscured  by  the  disk.  This  notion 
is  in  excellent  agreement  with  the  SIMTOI  model  and  observed 
photometry  (e.g.,  see  Figure  4).  A  bulk  comparison  of  the 
reconstructed  model  versus  real  data  implies  that  the  large 
concentration  of  flux  in  the  F-star’ s  northern  hemisphere  is 
probably  an  artifact,  whereas  the  small  crescent  of  the  F-star 
seen  in  the  east  is  real. 
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Table  9 

Best  fit  Values  and  Statistical  Information  for  Single  epoch  MultiNest  Minimizations  Involving  Three  Variants21  of  the  Tilted  Pascucci  Disk  Model  (Model  8)  The 

Posterior  Odds  Ratio  (A  log  R )  is  Relative  to  the  Pascucci  Zero  tilt  Disk  Model 


F  star 

Disk 

Statistical  Information 

#LDD(mas) 

«LDD 

t disk  (deg) 

^disk(deg) 

a 

rc(mas) 

p 

hc  (mas) 

up 

A  log  R 

X2r 

X2r(V2) 

xjffi  a) 

X2r(X34>) 

Model:  Pas 

cucci 

zero  tilt 

2009  11 

2.30 

0.63 

9.19 

1.87 

1.67 

0.023 

4708 

4.92 

4.23 

2.62 

8.34 

2009  12 

2.28 

1.00 

19.95 

3.71 

0.10 

0.118 

7046 

9.39 

4.00 

3.78 

25.02 

2010  02 

2.43 

0.98 

0.75 

0.173 

3046 

3.47 

2.83 

3.26 

4.66 

2010  08 

2.33 

0.74 

3.78 

0.011 

4349 

10.25 

4.00 

3.69 

26.18 

2010  09 

2.37 

0.72 

3.00 

0.021 

3336 

3.24 

1.98 

2.05 

6.48 

2010  10 

2.29 

0.66 

4.20 

0.007 

5569 

1.78 

1.33 

0.91 

3.49 

2010  11 

2.18 

0.10 

0.36 

0.093 

3106 

4.80 

2.91 

1.14 

11.29 

2010  12 

2.29 

0.41 

4.83 

0.005 

3726 

4.94 

2.92 

1.49 

11.84 

2011  01 

2.13 

0.18 

3.05 

0.015 

3030 

5.73 

3.31 

1.86 

13.71 

Average 

2.29 

0.60 

0 

0 

14.57 

2.79 

2.41 

0.052 

4215 

5.39 

3.06 

2.31 

12.34 

±  0.09 

± 

±  7.61 

±  1.30 

± 

dr  0.061 

dr 

0.32 

1.75 

1373 

Model:  Pas 

cucci 

fixed  tilt 

2009  11 

2.29 

0.62 

9.19 

1.79 

1.56 

0.026 

6329 

37 

4.91 

3.98 

2.64 

8.69 

2009  12 

2.28 

1.00 

19.94 

3.72 

0.10 

0.122 

6842 

119 

9.22 

3.93 

3.63 

24.66 

2010  02 

2.43 

0.99 

0.72 

0.179 

5 

296 

3.53 

2.92 

3.39 

4.59 

2010  08 

2.33 

0.74 

3.35 

0.014 

493 

2999 

10.65 

3.93 

3.82 

27.57 

2010  09 

2.39 

0.76 

2.81 

0.028 

129 

5802 

3.37 

2.02 

1.95 

6.98 

2010  10 

2.29 

0.66 

3.66 

0.011 

8 

1235 

1.70 

1.34 

0.88 

3.18 

2010  11 

2.18 

0.10 

0.37 

0.097 

29 

849 

4.69 

2.91 

1.15 

10.91 

2010  12 

2.29 

0.39 

4.80 

0.006 

56 

394 

5.22 

2.97 

1.56 

12.73 

2011  01 

2.13 

0.17 

2.27 

0.030 

81 

661 

5.50 

3.02 

1.68 

13.53 

Average 

2.29 

0.60 

1.3 

0.02 

14.56 

2.76 

2.18 

0.057 

4948 

55 

5.42 

3.00 

2.30 

12.54 

±  0.09 

± 

±  7.60 

±  1.36 

dr 

±  0.061 

± 

0.32 

1.62 

2245 

Model:  Pas 

cucci 

free  tilt 

2009  11 

2.30 

0.62 

1.102 

0.16 

8.93 

1.74 

1.55 

0.025 

6018 

349 

4.65 

4.37 

2.52 

7.25 

2009  12 

2.27 

1.00 

1.148 

0.15 

19.96 

3.73 

0.10 

0.120 

6417 

164 

9.14 

3.65 

3.29 

25.22 

2010  02 

2.43 

0.95 

1.088 

2.72 

1.10 

0.176 

46 

347 

3.06 

2.37 

2.63 

4.52 

2010  08 

2.33 

0.78 

2.330 

0.70 

4.75 

0.005 

1298 

1208 

16.17 

7.69 

7.96 

37.09 

2010  09 

2.39 

0.78 

1.153 

0.51 

2.70 

0.027 

422 

6353 

2.92 

1.71 

1.67 

6.11 

2010  10 

2.32 

0.71 

0.353 

0.29 

4.92 

0.003 

13 

1256 

4.43 

2.11 

1.19 

12.06 

2010  11 

2.24 

0.29 

6.879 

0.95 

0.24 

0.000 

691 

1511 

5.22 

2.93 

1.37 

12.52 

2010  12 

2.23 

0.38 

2.324 

0.96 

4.97 

0.000 

330 

780 

114.24 

22.38 

14.86 

370.25 

2011  01 

2.20 

0.23 

1.108 

0.87 

2.08 

0.037 

1391 

1971 

1.62 

1.25 

0.87 

3.01 

Average 

2.30 

0.64 

1.943  ± 

0.51 

14.45 

2.73 

2.49 

0.044 

4223 

522 

17.94 

5.39 

4.04 

53.12 

±  0.08 

db 

1.955b 

±  1.03 

±  7.80 

±  1.40 

dr 

dr  0.062 

dr 

0.28  1.97  1319 


Notes. 

a  The  variations  either  assume  the  disk  has  zero  tilt,  a  fixed  tilt,  or  a  per  epoch  tilt  with  respect  to  the  orbital  plane.  The  averaged  inclination  of  1  ?  33  ±  0?67  agrees 
well  with  the  multi  epoch  minimizations.  The  variations  in  scale  height  appear  to  be  real. 
b  Excluding  the  2010  11  result,  this  becomes  1.33  ±  0.67  in  agreement  with  our  fixed  tilt  model. 


A.7.  2011-09-18  and  Later 

All  of  the  post-eclipse  images  are  displayed  in  Figure  A52. 
Much  like  the  ingress  phase  images,  these  reconstructions  show 
that  the  F-star  has  no  egregious  asymmetries  that  cannot  be 
explained  by  UV  coverage  or  reconstruction  artifacts.  We  note 


that  two  of  the  post-eclipse  epochs  (2011-10-10  and  2011-11- 
03)  do  show  non-zero  closure  phase  and  different  locations  of 
the  first  visibility  null,  indicating  that  the  F-star  may  harbor 
spots  or  be  slightly  oblate.  Both  of  these  effects,  if  real,  are 
insignificant  compared  to  the  variations  that  the  disk  imparts 
upon  the  interferometric  data  during  the  eclipse.  We  will 
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Figure  5.  Best  fit  models  and  SQUEEZE  total  variation  (TV)  regularizer  reconstructions  for  each  CHAR  A  epoch.  All  images  are  displayed  with  North  up  and  east  to 
the  left.  The  0.5  mas  resolution  limit  of  MIRC  and  0.7  mas  limit  of  CLIMB  are  indicated  by  the  solid  white  circles  in  the  lower  left  of  the  first  image.  The  remaining 
white  circle  indicates  the  best  fit  SIM  TO  I  angular  diameter.  In  each  epoch  the  columns  contain  the  following  images  (left)  the  best  fit  model  from  SIMTOI,  (center) 
the  model  independent  reconstruction  of  the  interferometric  data,  and  (right)  a  reconstruction  of  synthetic  data  generated  from  the  model  as  described  in  Section  4.  For 
a  detailed  discussion  of  each  epoch  please  see  the  Appendix. 

(The  complete  figure  set  (7  images)  is  available.) 
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attempt  to  quantify  the  presence  of  spots  or  asymmetry  in 
greater  detail  in  a  future  publication. 

APPENDIX  B 

CREATION  OF  THE  SILHOUETTE  IMAGE 

The  silhouette  image  is  created  by  aligning  and  combining 
the  individual  SQUEEZE  images.  Because  SQUEEZE  images 
are  flux  centered  and  frequently  contain  stray  flux,  the  image 
was  generated  manually  using  the  following  process.  SIMTOI 
was  used  to  generate  a  series  of  768  x768  pixel  images 
consisting  of  a  uniform  disk  F-star  aligned  according  to  the 
orbital  elements  specified  in  this  work.  These  images,  along 
with  the  original  SQUEEZE  reconstructions,  were  imported 
into  GIMP  (a  popular  image  editing  program)  as  layers.  During 
the  import  process,  each  FITS  image  was  converted  to  256 
levels  of  gray.  Individual  SQUEEZE  layers  were  aligned  with 
the  corresponding  SIMTOI  prediction  with  the  aid  of  a  ranged 
color  selection  tool  used  to  highlight  the  average  background 
plus  6%.  The  layers  were  combined  using  the  “screen”  layer 
mode  which  computes  the  output  value,  C,  from  overlays  A  and 
B  according  to  the  following  equation: 

C  =  255  -  [(255  -  A)  •  (255  -  B)]/ 255. 

Because  each  image  is  of  unit  flux,  this  scaling  approach 
preserves  the  appearance  of  the  original  images.  After  this,  a 
black  mask  was  created  that  decreased  the  appearance  of  stray 
flux  in  the  individual  images.  This  last  step  was  subjective  and 
therefore  used  sparingly.  Finally,  the  image  was  gamma- 
corrected  7  =—1.5  in  order  to  make  the  appearance  of  the 
southern  limb  of  the  F-star  more  obvious  to  the  reader. 
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